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SECTION  I 


INTRODUCTION  AND  SUMMARY 


An  essential  requirement  of  military  STOL  tactiaal  transports 
planned  for  the  1980  time  period  will  be  to  operate  from  airfields 
of  2500  feet  or  less.  These  aircraft  will  use  thrust  reversers 
as  primary  braking  devices  throughout  the  landing  ground  roll. 

Also,  some  STOL  concepts  will  use  thrust  vectoring  systems  to 
help  control  the  flightpatn  of  the  airplane  and  reduce  take¬ 
off  and  landing  speeds.  Consequently,  emphasis  must  be  placed  on 
designing  efficient  and  reliable  thrust  reverser/vectoring  systems 
to  achieve  the  field  length  objective. 

Commercial  jet  aircraft  have  used  thrust  reversers  as  secondary 
braking  devices  since  the  beginning  of  their  operation.  However, 
the  complex  problems  caused  by  the  interactions  between  reverser 
exhaust  and  aircraft  flowfields  have  limited  their  usefulness. 

These  problems  include  exhaust  gas  recirculation  which  can  lead 
to  engine  surge,  impingement  of  exhaust  gases  on  the  ground  or 
adjacent  aircraft  surfaces,  and  engine  mass  flow  matching.  Also, 
the  reverser  flow  can  cause  blanking  out  of  aerodynamic  control 
surfaces  leading  to  a  loss  in  aircraft  directional  stability 
and  control,  buoyancy  effects  that  decrease  the  efficiency  of 
the  ground  braking  systems,  and  changes  in  airplane  drag.  All 
of  these  problems  have  been  experienced  during  the  development 
of  existing  commercial  aircraft.  However,  the  availability  of 
long  runways  has  made  it  unnecessary  to  completely  resolve  the 
interactions  between  the  reverser  and  aircraft  flowfields. 

To  avoid  the  limitations  of  existing  systems  on  future  STOL  aircraft, 
attention  must  be  given  to  the  following  technical  areas: 

o  TR/TV  performance 

o  Exhaus4-  gas  flowfield 

o  Aerodynamic  interference 

o  Engine  operation 

o  TR/TV  system  design  including  weights  and  structures 

The  above  considerations  have  significant  influence  on  nacelle 
placement,  thrust  reverser  and  vectoring  system  geometry,  and 
operating  envelope. 

The  Boeing  Company,  with  subcontract  support  from  Pratt  &  Whitney 
Aircraft  conducted  an  18  month  research  program  to  study  the  above 
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technical  areas.  The  program  was  administered  by  the  Air  Force 
Aero  Propulsion  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio. 
Program  objectives  are: 

1.  To  develop  methods  to  predict  thrust  reverser  and  thrust 
vectoring  system  performance. 

2.  To  establish  design  criteria  for  high  efficiency,  lightweight 
thrust  reversers  or  thrust  vectoring  systems  for  STOL  aircraft. 

The  program  has  three  parts: 

Part  IA — Data  Review  and  Analysis 

Part  IB — Design 

Part  IC — Model  Testing 

Part  IA  consists  of  three  tasks: 

Task  1.1  -  Review  and  correlate  TR/TV  data 

Task  1.2  -  Construct  computerized  analytical  models 

Task  1.3  -  Plan  and  conduct  supplemental  tests 

During  Task  1.1,  existing  data  were  reviewed  for  possible  appli¬ 
cation  to  computer  programs  for  thrust  reverser  and  vectoring  systems. 
Literature  searches  of  DDC,  NASA,  and  Boeing  files  resulted  in 
approximately  160  references  applicable  to  TR/TV  systems.  The  re¬ 
sults  of  the  data  review  were  used  to  develop  data  correlations  for 

o  Cruise  nozzles 

o  Thrust  reversers 

o  Thrust  vectoring  nozzles 

o  Cascade  lattices  applicable  to  cascade  TR  and  TV  nozzles 

Computer  programs  for  predicting  TR  and  TV  nozzle  performance  and 
evaluating  TR  and  TV  influence  on  the  total  airplane  system  were 
developed  during  Task  1.2.  Three  programs  were  developed: 

1)  Jet  Trajectory  and  Spreading  Program  —  to  predict  the  shape 

and  trajectory  of  the  thrust  reverser  or  vectoring  nozzle  exhaust 
plume. 

2)  Reingestion  Prediction  Program  —  to  predict  the  on-set 
of  reingestion  for  arbitrary  thrust  reverser  and  air¬ 
plane  configurations  as  a  function  of  geometry  and  flow 
conditions. 

3)  TR  and  TV  System  Performance  Program  —  consisting  of  four 
modules  to  predict: 
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o  TR  and  TV  Internal  Performance 
o  Aerodynamic  Performance 
o  Reingestion 
o  Engine  Stability  Margin 

The  TR  and  TV  Internal  Performance  Module  was  assembled  using 
the  data  correlations  developed  during  Task  1.1.  The  Engine 
Stability  Margin  Module  was  developed  by  Pratt  &  Whitney 
Aircraft.  Available  reingestion  and  aerodynamic  data  for  existing 
conventional  take-off  and  landing  aircraft  were  not  applicable 
to  the  Reingestion  and  Aerodynamic  Performance  nodules.  A  low 
speed  wind  tunnel  test  of  STOL  airplane  TR  and  TV  configurations 
is  required  to  obtain  the  necessary  data.  However,  logic  was 
provided  to  allow  easy  incorporation  of  data  into  the  modules 
as  data  become  available. 

Task  1.3  consisted  of  planning  and  conducting  supplemental  static 
test  to  fill  data  voids  discovered  in  the  open  literature.  Tests 
were  conducted  to  determine: 

1)  Multibearing  vectoring  nozzle  performance  as  a  function  of 
parametric  geometry  variations 

2)  Blocker  door  geometry  effects  on  the  performance  of  annular 
blocker/cascade  thrust  reversers 

The  results  were  incorporated  into  the  Internal  Performance  Module 
of  the  TR  and  TV  System  Performance  Program. 

Detailed  descriptions  of  the  results  of  Part  IA  are  provided  in 
the  following  sections.  Detailed  results  of  Parts  IB  and  IC 
are  provided  in  Volume  II. 
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SECTION  II 


PART  IA  -  DATA  REVIEW  AND  ANALYSIS 


2.1  Task  1.1  Review  and  Correlate  TR/TV  Data 


The  objective  of  Task  1.1  was  to  review  the  existing  literature 
for  data  pertinent  to  TR/TV  systems  and  to  correlate  the  data 
as  a  function  of  fundamental  geometric  and  aerodynamic  para¬ 
meters.  The  data  correlations  are  used  to  predict  TR  and  TV 
nozzle  performance.  The  first  step  of  this  task,  a  litera¬ 
ture  search,  was  made  prior  to  start  of  the  study  contract.  The 
literature  search  included  the  following  sources  of  information. 


1) 

Defense  Documentation  Center  computerized 
literature  search.  This  search  was  updated 
during  Task  1.1 

2) 

Computerized  literature  search  of  NASA  reports 

3) 

Boeing  documents,  STAR  and  TAB  abstracts, 
technical  journals  (search  performed  by 

Boeing  library  personnel) 

4) 

Foreign  literature  available  through  services  of 
the  Boeing  International  Corporation 

5) 

Patent  search 

Approximately  160  reports  found  during  the  literature  search 
were  obtained  and  reviewed.  Boeing  was  assisted  in  the  liters 
ture  review  by  the  subcontractor,  Pratt  &  Whitney  Aircraft, 

who  provided  references  and  abstracts.  Results  were  publish¬ 
ed  as  an  Air  Force  Technical  Report  (Ref.  1).  The  data  review 
document  contains  three  types  of  information. 

Bibliography  Summary  Chart 

A  typical  bibliography  summary  chart  for  thrust  reverser 
systems  is  shown  in  Table  I.  This  chart  summarizes  information 
extracted  from  reports  and  cross  references  them  by  subject 
matter. 

Report  Abstracts 

Abstracts  were  written  for  all  reports  reviewed.  Abstracts 
describe  contents  of  the  report  and  provide  an  objective 
assessment  of  applicability  and  usefulness  of  the  contents. 
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Data  Review  Charts 


A  sample  data  review  chart  is  shown  in  ^able  II.  This  chart 
contains  a  sketch  of  the  TR  or  TV  configuration ,  summarizes 
the  range  of  test  variables ,  and  lists  the  type  of  data  con¬ 
tained  in  the  report.  Comments  are  included  conoerning  the  model, 
data  quality^ or  usefulness  of  the  data. 

The  second  major  objective  of  Task  1.1  was  to  formulate  data 
correlations  for  thrust  reverser  and  thrust  vectoring  systems 
as  functions  of  fundamental  geometric  and  aerodynamic  parameters. 
Correlations  were  developed  for  the  following  types  of  TR/TV 
nozzles: 

1)  Cruise  nozzles 

a)  Conical 

b)  Annular 

c)  Noncircular 

2)  Thrust  reversers 

a)  Target  (clamshell  and  annular) 

b)  Blocker  deflector  and  blocker  cascade 

3)  Thrust  vectoring  nozzles 

a)  Single  bearing 

b)  Three  bearing 

c)  Spherical  eyeball 

d)  Lobs ter tail 

e)  External  deflector 

4)  Cascade  lattices  applicable  to  cascade  TR  and  TV  nozzles 

Data  correlations  for  the  above  TR/TV  nozzles  are  described 
in  the  following  paragraphs. 

2.1.1  Cruise  Nozzle  Data  Correlations 

This  section  describes  methods  used  to  predict  cruise  nozzle 
performance  in  terms  of  velocity  and  discharge  coefficients. 
Extensive  parametric  data  for  convergent  conical  nozzles  are 
presented  covering  a  wide  range  of  wall  angles  and  diameter 
ratios.  Analytical  results  are  also  presented  tnd  compared 
to  experimental  data. 

Conical  Nozzle  Velocity  Coefficient 

A  sketch  of  a  conical  nozzle  is  shown  on  page  8. 
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Table  II:  TYPICA  L  THRUST  RE  VERSE R  DA  TA 

REVIEW  CHART  FOR  CASCADE  THRUST 
REVERSERS 


Results  from  a  parametric  test  of  16  convergent  conical 
nozzles  (Ref.  2)  are  shown  in  Figures  1  to  4.  The  data  show 
small  but  definite  effects  due  to  wall  angle  0(  and  diameter 
ratio  D.,/D2.  These  data  are  used  in  the  Internal  Performance 
Module  *of  the  TR  and  TV  System  Performance  Program.  Pen¬ 
alties  due  to  skin  friction  and  underexpansion  losses  are 
charged  separately,  as  shown  in  Figure  1.  If  the  nozzle  exit 
is  offset  from  the  nacelle  centerline,  which  is  typical  of 
single  bearing  vectoring  nozzle  designs,  then  the  Acv  penalty 
shown  in  Figure  5  is  charged. 

During  preliminary  design  studies,  the  simplest  method  of  pre¬ 
dicting  nozzle  performance,  termed  a  "Level  1"  prediction,  is 
adequate  because  the  nozzle  geometry  is  not  well  defined.  The 
nozzle  is  charged  with  Ac„  =  0.005  for  skin  friction  losses 
and  offset  losses  from  Figure  5.  In  addition,  underexpansion 
losses  are  charged  for  supercritical  pressure  ratios  as  shown 
in  Figure  6. 


Conical  Nozzle  Discharge  Coefficient 

Experimental  results  for  convergent  nozzle  discharge  co¬ 
efficients  (Ref.  2)  are  shown  in  Figure  7.  Note  that  wall 
angle,  diameter  ratio,  and  pressure  ratio  have  significant 
effects  on  discharge  coefficient.  Experimentally  determined 
choking  pressure  ratios  are  indicated  on  the  curves  by  tic 
marks.  A  shaded  band  has  been  drawn  bracketing  the  experimental 
points.  The  choked  discharge  coefficient  levels  are  summariz¬ 
ed  in  Figure  8  as  a  function  of  wall  angle  and  diameter 
ratio. 


Annular  Nozzle  Velocity  Coefficient 
A  sketch  of  an  annular  nozzle  is  shown  below. 
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Figure  5:  EFFECT  OF  NOZZLE  OFFSET  ON  VELOCITY 
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13 


A0  'iNilDHdaOO  AXI3013A 


14 


Figure  6:  THEORETICAL  VELOCITY  COEFFICIENT 

UNDEREXPANSION  LOSSES  FOR  CONVERGENT 
NOZZLES  AT  SUPERCRITICAL  PRESSURE  RATIOS 


1.0 


3.0 


2.0 

NOZZLE  PRESSURE  RATIO.  PTN/P«> 

Figure  7:  EXPERIMENTAL  DISCHARGE  COEFFICIENT 

CURVES  FOR  CONVERGENT  CONICAL  NOZZLES 


15 


FfanB:  PARAMETRIC  STUDY  OF  CONICAL  CONVERGENT 
NOZZLES  -  CD  CHOKE 


Test  results  are  used  to  predict  the  Acv  penalty  of  annular 
nozzles  relative  to  conical  nozzles.  Results  from  Boeing 
tests  of  JT3D  and  C-5A  fan  nozzles  are  shown  in  Figure  9.  The 
conical  nozzle  Cv  data  are  charged  with  this  additional  increment 
to  predict  annular  nozzle  performance. 


Annular  Nozzle  Discharge  Coefficient 


Review  of  the  literature  did  not  produce  a  usable  set  of  data 
for  annular  nozzle  discharge  coefficients.  However,  a  combin¬ 
ation  of  conical  nozzle  data  and  theory  was  developed  that  is 
adequate.  Boundary  layer  theory  is  used  to  determine  discharge 
coefficient  assuming  only  viscous  losses,  i.e.,  no  vena  con- 
tracta  or  three-dimensional  flow  effects.  For  a  conical  nozzle, 

.  l-lf 


o  Conic*  I  nooclc 

where  R  is  the  nozzle  exit  radius  and  8*  is  the  boundary 
layer  displacement  thickness  at  the  nozzle  exit.  For  an 
annular  nozzle, 

c  . 

0 annular-  n*«»lc 


(1) 


(2) 


where  R.and  R  are  the  annular  nozzle  inner  and  outer  radii. 
Assuming  that°the  conical  and  annular  nozzles  have  the  same 
length,  area,  exit  Mach  number,  and  Reynolds  number 
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« 

Conical 


°  annular 
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Eliminating  8*  between  Equations  1  and  2  and  solving  for 
Sannular  «ives 


C 


O  annular 


R.-R  a 


0-c 


Ocom'cal' 


(4) 


The  term  cDconicai  is  obtained  from  the  conical  nozzle  data 

at  choked  flow  conditions,  Figure  8.  Equation  4  is  used  to 
calcdate  C’D  annu^ar  at  choked  conditions. 


In  order  to  predict  C_  at  subcritical  pressure  ratios,  the 
analytical  method  of  Bragg  (Ref.  3)  is  used.  Bragg's  analytical 
discharge  coefficient  curves  are  shown  in  Figure  10.  Bragg's 
results  are  erroneous  for  pressure  ratios  greater  than  choke. 
Consequently,  the  curves  were  corrected  by  making  C  constant 
above  the  choked  pressure  ratio.  The  choking  line  was  taken 
from  Figure  7  .  A  good  check  of  the  accuracy  of  Bragg's  theory 
is  made  by  overlaying  Figures  7  and  10. 

Bragg’s  theory  does  not  permit  the  discharge  coefficient  to  be 
calculated  for  a  particular  nozzle  geometry.  It  predicts  how 
the  discharge  coefficient  of  a  nozzle,  known  under  choked  flow 
conditions  will  vary  for  subcritical  pressure  ratios.  To  use 
Bragg's  method,  the  choked  discharge  coefficient  is  located  on 
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Figure  9:  EFFECT  OF  NOZZLE  OFFSET  ON  ANNULAR  NOZZLE 
VELOCITY  COEFFICIENT  LOSSES 
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DISCHARGE  COEFFICIENT.  C 


the  right  side  of  Figure  10.  The  curve  is  followed  back  to  the 
left  by  fairing  in  the  best  fit  at  the  choke  line. 


Noncircular  Nozzle  Velocity  Coefficient 

Noncircular  nozzles  have  rectangular,  O-shaped,  or  elliptical 
shaped  cross  sections  as  shown  in  the  sketches  below. 


( - 1 

A 

^  1 

_ / 

The  prediction  method  for  noncircular  nozzles  employs  an 
equivalent  hydraulic  diameter  correlation  developed  to  predict 
the  maximum  velocity  coefficient  of  suppressor  nozzles,  as 
shown  in  Figure  11  from  Ref.  4.  If  the  losses  are  assumed  to  be 
caused  predominantly  by  skin  friction  losses  in  the  boundary 
layer,  the  maximum  velocity  coefficient  of  circular  convergent 
nozzle  is  given  by 


C.ire  uUr 


4-0 

D“ 


(5) 


where  0  is  the  diameter  and  0  is  boundary  layer  momentum 
thickness  at  the  nozzle  exit.  The  maximum  velocity  coefficient 
occurs  at  a  pressure  ratio  such  that  the  convergent  nozzle  is 
choked  but  does  not  exhibit  any  underexpansion  loss.  This 
typically  occurs  at  a  nozzle  pressure  ratio  of  about  2.2.  By 
dimensional  analysis,  it  can  be  shown  that  for  noncircular 
nozzles, 


VMAk 


noncircuiar 


(6) 


where  C  is  the  maximum  velocity  coefficient  of  a  standard 
convergent  nozzle,  and  D^e  is  the  equivalent  hydraulic  diameter 
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Cy  max. MAXIMUM  VELOCITY  COEFFICIENT 
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Figure  1 1:  MAXIMUM  VELOCITY  COEFFICIENT  CORRELA  TION 
FOR  SUPPRESSOR  NOZZLES 
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The  percent  loss  in  velocity  coefficient  relative  to  the 
standard  nozzle  is  given  by 


A Cv  _  I  _  jocular  ^gj 

CvS  Cyg 

To  apply  the  equivalent  hydraulic  diameter  correlation  to  a 
noncircular  nozzle,  the  parametric  velocity  coefficient  data 
shown  in  Figures  1  to  4  are  used  to  obtain  the  nozzle  skin 
friction  loss  Ac„.  The  increase  in  skin  friction  loss  due  to 
increased  wetted  area  of  the  noncircular  nozzle  is  given  by 

A  r  _  (i  (CvmAn)  \  (  i  Ar  \ 

arm  ~  \  ~  Cvs  /  V  "  frfctiW  (9) 


Noncircular  Nozzle  Discharge  Coefficient 

Boundary  layer  theory  is  employed  to  predict  discharge  co¬ 
efficient  penalties.  Assuming  only  skin  friction  losses. 


C 


v  «kl  A  frictiori 
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Combining  Equations  10  and  11  gives 

^  ^0  Skin  fricWo*  *  ^  ^  ^ V  skin  fWcti**  (12) 

where  H  «  S*/0  is  boundary  layer  shape  factor.  For  a  turbulent 
flatplate  boundary  layer  with  M,,  ■  1,  a  good  value  for  H  *  1.7. 
The  discharge  coefficient  penalty  due  to  wetted  area  increase 
is  given  by 


AC 
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The  maximum  discharge  coeeficient  is  given  by 


(13) 


*  ^Oc»*ka  ”  ACe  ai>«  (14) 

where  CD  c^0^e  is  obtained  from  the  conical  nozzle  choked  dis¬ 
charge  coefficient  data,  Figure  8.  Discharge  coefficient  is  pre 
dieted  at  lower  pressure  ratios  using  Bragg's  theoretical 
curves  with  CDjnax  from  Equation  14  as  the  choked  discharge  co¬ 
efficient  value. 
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2.1.2  Thrust  Reverser  Data  Correlations 

This  section  describes  data  correlations  developed  for  several 
types  of  thrust  reverser  systems.  Boeing  was  assisted  in 
|  formulating  the  data  correlations  by  Pratt  &  Whitney  Aircraft 

who  provided  data  for  annular  target  and  blow-in  door 
ejector  thrust  reverters. 
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Clamshell  Target  Thrust  Reverser 


A  clamshell  target  thrust  reverser  is  shown  in  Figure  12.  Clam¬ 
shell  target  reversers  are  used  on  the  737,  DC-9,  and  C-141  air¬ 
planes.  An  experimental  study  of  geometric  variables  influenc¬ 
ing  the  static  performance  of  clamshell  target  reversers  is 
presented  in  Ref.  5.  The  geometric  variables  are  defined  in 
Figure  13.  Table  III  summarizes  the  range,  nominal  value,  and 
the  figure  number  showing  the  effect  of  each  geometric  variable 
on  corrected  static  reverser  efficiency. 


Table  III:  SUMMARY  OF  CLAMSHELL  TARGET  THRUST  REVERSER 
I  GEOMETRIC  VARIABLES 


In  addition  to  the  geometric  variables,  nozzle  pressure  ratio  was 
varied  from  1.2  to  1.98.  Data  were  taken  for  corrected  static 


reverser  efficiency,  and  airflow  match,  §.  In  functional 
form,  Rc 


c  (  Ptm  *  l  in. 
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TARGET  DOORS 
(DEPLOYED  POSITION) 


Ffan  12:  CLAMSHELL  TARGET  THRUST  REVERSER 
SCHEMATIC 


I 


> 


Figure  13:  CLAMSHELL  TARGET  THRUST  REVERSER 
GEOMETRIC  VARIABLES 
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CORRECTED  STATIC  REVERSER  EFFICIENCY.?/  RC  «  (Fx/W)REV/<F/W»F, 


DOOR  LENGTH  RATIO 
AVERAGE  LIP  HEIGHT  = 
SWEEP  ANGLE  =■  X  =  10° 
ARC  ANGLE  =  0  =  140° 
CONE  ANGLE  =  a  =  10° 
BEVEL  ANGLE  =  0  =  0° 
NO  SIDE  FENCES 


L/D  =  0.95 
LH/D  =  0.106 


0.B  0.9 

SETBACK  RATIO,  X/D 

Figun  14:  EFFECT  OF  SETBACK  RATIO  ON  STATIC 
REVERSER  EFFICIENCY,  CLAMSHELL 
TARGET  THRUST  REVERSER 
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CORRECTED  STATIC  REVERSER  EFFICIENCY,  i?RC=  (Fx^REV^^FWD 


SETBACK  RATIO  »  X/D  =  0.S4 
AVERAGE  LIP  HEIGHT  -  LH/D  =*  0.106 


SWEEP  ANGIE  = 

X  = 

10° 
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0  » 

140° 

CONE  ANGLE  - 

a  * 

10° 

BEVEL  ANGLE  » 

0  - 

0° 

NO  SIDE  FENCES 

Figure  IS:  EFFECT  OF  DOOR  LENGTH  ON  STA  TIC 
REVERSER  EFFICIENCY.  CLAMSHELL 
TARGET  THRUST  REVERSER 
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CORRECTED  STATIC  REVERSER  EFFICIENCY,  tjRC  -  <FxAWREV/(F/W)F 
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CORRECTED  STATIC  REVERSER  EFFICIENCY,  r?RC  =  (Fx/W)rev/F/W)FWi:) 


SETBACK  RATIO  =  X/D  =  0.94 
DOOR  LENGTH  RATIO  *  L/D  *1.0 
AVERAGE  LIP  HEIGHT  =  LH/D*  0.106 
ARC  ANGLE  =  6  =  140° 

CONE  ANGLE  =  a  =  10° 

BEVEL  ANGLE  =0  =  0° 

NO  SIDE  FENCES 


Figure  17:  EFFECT  OF  SWEEP  ANGLE  ON  STATIC  REVERSER  EFFICIENCY, 
CLAMSHELL  TARGET  THRUST  REVERSER 
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CORRECTED  STATIC  REVERSER  EFFICIENCY.  t?RC  -  IF 


SETBACK  RATIO  -  X/D  -  0J4 
DOOR  LENGTH  RATIO  -  Up  -  lAi 
AVERAGE  UP  HEIGHT  -  LH/D  •  0.106 
SWEEP  ANGLE  -  X  -  10*' 

CONE  ANGLE  -  a  -  10° 

BEVEL  ANGLE  -0*0 P 
NO  SIDE  FENCES 


Figun  18:  EFFECT  OF  ARC  ANGLE  ON  STA  TIC  REVERSER 
EFFICIENCY,  CLAMSHELL  TARGET  THRUST 
REVERSER  « 


SETBACK  RATIO  -  X/D  -  034 
DOOR  LENGTH  RATIO  -  L/D  -  1.0 
AVERAGE  LIP  HEIGHT  «  IB/D  *  0.106 
SWEEP  ANGLE  «  X  -  10° 

ARC  ANGLE  *  6  -  140° 

BEVEL  ANGLE  ■=  0  =  0° 


-  »  - - f- 

0  5  10 

TARGET  DOOR  CONE  ANGLE,  a,  DEGREES 


Figure  19:  EFFECT  OF  CONE  ANGLE  ON  STATIC  THRUST 
REVERSER  EFFICIENCY,  CLAMSHELL  TARGET 
THRUST  REVERSER 
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CORRECTED  STATIC  REVERSER  EFFICIENCY.  r?RC  =  <fyW)REV/<F/W>F 


SETBACK  RATIO  »  X/D  *  0.04 
DOOR  LENGTH  RATIO  «  JL/D  *  1.0 
AVERAGE  UP  HEIGHT  -  LH/D  -  0.106 
SWEEP  ANGLE  -  X  -  10° 

ARC  ANGLE  -  6  -  140° 

CONE  ANGLE  -  a  -  10° 

NO  SIDE  FENCES 


Figurt  20:  EFFECT  OF  BEVEL  ANGLE  ON  STA  TIC 
REVERSER  EFFICIENCY,  CLAMSHELL 
TARGET  THRUST  REVERSER 
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Nozzle  pressure  ratio  and  setback  distance  X/D  were  the  only 
variables  that  significantly  affected  airflow  match. 


(16) 


Effects  of  setback  distance  on  airflow  match  are  shown  in  Figure 
21.  The  baseline  corrected  reverser  efficiency  and  airflow  match 
curves  are  shown  in  Figure  22.  The  procedure  used  to  predict  YRc 
and  i  for  an  arbitrary  clamshell  geometry  is  to  assume  that 
effects  of  the  variables  in  Equation  15  are  independent.  This 
assumption  is  necessary  because  the  test  did  not  include  all 
possible  combinations  of  parameters.  The  equations  for  corrected 
reverser  efficiency  and  airflow  match  then  become: 


=  Rc)gA5£UN£^  (^*c)>f/0  *  Yd  +  (^*c)lh/d 

+  (A7J,C),  4  (A^  +(AW  +(A^c)a  (17) 

(A^v/o  U8) 

where  the  incremental  terms  are  found  by  subtracting  the  baseline 
values  from  Figures  14  to  21,  e.g. 

(A^Rt)x/b  ‘  (^Ox/o  *  (^JsAstUNt  (19) 

'  W  FIG.  14  flQ.Zl 

Annular  Target  Thrust  Reverser 

A  sketch  of  an  annular  target  thrust  reverser  is  shown  below. 


Several  excellent  data  sources  exist  for  annular  reversers 
(Ref.  6,  7  &  3).  Acorrelation  for  corrected  reverser 

efficiency  that  summarizes  the  data  sources  is  shown  in  Figure 
23.  The  correlating  parameter  for  is  blockage  angle  minus 

door  angle,  P-9  ,  as  identified  in  Figure  23.  As  P-9  increases, 
^Rc  aPProac^es  the  limiting  value  of  cos  6  .  The  parameter 
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0.9 


04 


SETBACK  RATIO,  XID 


1 


Figun21:  EFFECT  OF  SETBACK  RATIO  ON  AIRFLOW 

MATCH,  CLAMSHELL  TARGET  THRUST  REVERSER 
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CORRECTED  REVERSE  THRUST/CORRECTED  AIRFLOW  RATIO. 


Figure  21  EFFECT VF  BLOCKAGE  AND  DOOR  ANGLE 
ON  CORRECTED  REVERSER  EFFICIENCY 
ANNULAR  TARGET  THRUST  REVERSER 
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pf-e  was  found  to  improve  the  ability  to  correlate  annular 
TR  data. 

The'  airflow  match  data  are  correlated  as  a  function  of  throat 
gap  over  annulus  height  in  Figure  24.  The  Boeing  707  and  Pratt 
&  Whitney  Aircraft  data  are  for  a  pos.texit  annular  target 
reverser.  Because  the  controlling  area  for  a  postexit  reverser 
is  the  cruise  nozzle  arei, the  limiting  value  for  airflow  match 
i  equals  1.0.  On  the  other  hand,  the  Boeing  C-5A  proposal  data 
are  for  a  pre-exit  annular  blocker  deflector  reverser.  The 
controlling  area  in  the  TR  mode  is  larger  than  the  cruise  nozzle 
are  so  i  is  greater lhan  1.0.  The  equation  for  airflow  match  is 

r  _  w r  _  Cpi»  A* 

wp  Cof  Ap  (20) 

where  w  is  mass  flow  rate,  CD  is  discharge  coefficient,  A  is 
exit  area  and  subscripts  r  and  f  refer  to  reverser  and  forward 
thrust  modes,  respectively.  Equation  (20)  was  used  to  calculate 
discharge  coefficient  ratio  from  the  Boeing  C-5A  airflow  match 
data.  The  resulting  data  correlation  curve  for  discharge  co¬ 
efficient  ratio  is  shown  in  Figure  25. 


Blocker  Deflector  Thrust  Reverser 

Data  correlations  were  developed  for  two  types  of  blocker  de¬ 
flector  reversers,  a  blocker  with  cascade  deflectors,  and  a 
blocker  reverser  with  deflector  doors.  Cascade  blocker  deflector 
thrust  reversers  are  used  on  the  707,  747,  DC-10,  and  L-1011 
airplanes.  Blocker  reversers  with  deflector  doors  are  used  on 
the  727  and  F-11A  airplanes. 

The  cascade  blocker  deflector  shown  in  Figure  109  on  page  163 
was  tested  during  the  Task  1.3  supplemental  static  test  and  the 
results  are  dit cussed  in  Section  2.3.  Baseline  performance  of 
the  blocker  reverser  with  deflector  doors  is  presented  in  Figure 
26  from  Ref.  9.  The  effects  of  blocker  door  cone  angle  on 
reverser  efficiency  and  airflow  match  are  displayed  in  Figure  27. 

2.1.3  Thrust  Vectoring  Nozzle  Data  Correlations 

This  section  describes  data  correlations  developed  for  several 
types  of  thrust  vectoring  nozzles.  Boeing  was  assisted  in  this 
ask  by  Pratt  &  Whitney  Aircraft,  who  provided  data  for  single- 
faring  and  spherical  eyeball  vectoring  nozzles. 
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THROAT  GAP/ANNULUS  HEIGHT,  C/H 

Figure  24:  EFFECT  OF  DOOR  SETBACK  ON  AIRFLOW 

MATCH  ANNULAR  TARGET  THRUST  REVERSER 
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THROAT  GAP/ANNULUS  HEIGHT,  C/H 
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TARGET  THRUST  REVERSE R 


Figun  26:  BASELINE  PERFORMANCE  FOR  BLOCKER-DEFLECTOR 
THRUST  REVERSER 
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NOZZLE  PRESSURE 


AIRFLOW  MATCH.  4»  =  — STATIC  REVERSER  EFFICIENCY 


Single-Bearing  Nozzle 

Sketches  of  single-bearing  nozzles  are  shown  below 


Considerable  geometric  flexibility  is  possible  by  varying  the 
nozzle  offset ,  bearing  plane  angle ,  and  bearing  duct  angle. 
Single-bearing  nozzles  frequently  are  designed  in  symmetric 
pairs  to  minimize  assymetric  engine  side  loads.  Theoretical 
thrust  components  for  a  single-bearing  nozzle  or  half  of  a 
dual  swiveling  single-bearing  nozzle  are  given  by 


Fx/Fr*  I  -  Sin ^43  (l-COS$ 


(21) 


Fy/Fr  *  SINK  SWtf  COS/3  (l-COS^)  -  C05K  SINtfSIN^ 


(22) 


F*/Fr  *  -COSH  SIN£  COSfi  ( I-  COS^  -  SINK  S INfi  S 
where  the  resultant  force  Fr  is  given  by 
fr,  (fxl  4  F/  4F4l)/l 


(23) 


(24) 


The  (x,  y ,  z)  coordinate  system  and  angles  0(,  $,  and  f  are  defined 
in  Figure  28.  Theoretical  thrust  components  for  a  nozzle  with 
bearing  plane  and  duct  angles  of  45  and  58  degrees  are  given  in 
Figure  29.  The  theoretical  flow  turning  angle  is  given  by 


A 


- 1 

COS 


(25) 
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DUCT  CENTERLINE 


NOZZLE  CENTERLINE 


PLANE  ANGLE 
DUCT  ANGLE 
ROTATION  ANGLE 


Figure  28:  SINGLE  BEARING  NOZZLE  NOMENCLA  TURE 
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Vector  angle  is  given  by  the  following  equation 


CT 


(26) 


Several  data  sources  exist  for  single-bearing  nozzle  performance 
(Ref.  8,  lo,and  n) .  As  shown  in  Figure  30,  the  various  data 
sources  agree  fairly  well  considering  the  geometric  differences 
between  test  configurations.  Consequently,  data  supplied  by 
Pratt  &  Whitney  Aircraft  was  selected  for  use  in  the  Internal 
Performance  Module. 
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l'hree-Bearing  Nozzle 


Data  measured  during  the  Task  1.3  supplemental  static  tests  are 
used  to  predict  three-bearing  nozzle  performance.  The  data 
are  discussed  in  Section  2.3. 


I  Spherical  Eyeball  Nozzle 

A  sketch  of  a  spherical  eyeball  nozzle  is  shown  below. 


Velocity  and  discharge  coefficient  performance  is  shown  in 
Figure  31  from  Ref.  8  and  11. 

Lobstertail  Nozzle 

A  lobstertail  (also  known  ?.s  aft  hood  deflector)  sketch  is  shown 
below. 
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VELOCITY  COEFFICIENT.  Cw  DISCHARGE 


Figure  30:  SINGLE  BEARING  NOZZLE  PERFORMANCE 
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VELOCITY  COEFFICIENT.  Cy  DISCHARGE  COEFFICIENT,  C 


DISCHARGE  COEFFICIENT.  CD,  OR  VELOCITY  COEFFICIENT.  C 


Velocity  and  discharge  coefficient  performance  for  95-degree 
vector  angle  is  shown  in  Figure  32  from  Ref.  8.  Effects  of 
vector  angle  and  nozzle  pressure  ratio  on  Cy  and  CD  are  given  in 
Figure  33.  The  cruise  nozzle  Cv  and  CD  arevestimated  to  be 
0.995  and  0.98,  respectively  atvchokeduconditions .  Bragg's 
theoretical  discharge  coefficent  curves  are  used  for  CD  at 
intermediate  vector  angeles. 


External  Deflector  Nozzle 


External  deflector  nozzles  employ  a  flat  plate  or  curved  surface 
downstream  of  the  nozzle  exit  to  deflect  the  flow,  as  depicted 


"hinged" 


Variables  include  nozzle  pressure  ratio,  deflection  angle  6, 
setback  ratio  X/D,  and  door  length  ratio  L/D.  The  effect  of 
deflection  angle  on  Cv  and  choked  CQ  for  flat  plate  and  curved 
deflectors  is  shown  in  Figure  34  for  nozzle  pressure  ratio 

0.  Choked  discharge  coefficient  is  used  with  Bragg's 
curves  to  give  C_  at  subcritical  pressure  ratios. 

The  effect  of  pressure  ratiouon  Cy  at  70°  deflection  angle  is 
shown  in  Figure  35. 

Existing  data  were  found  inadequate  to  predict  effects  of 
setback  and  door  length  ratio.  Theoretical  results  showing  the 
effect  of  door  length  ratio  on  deflection  angle  are  shown  in 
Figure  36  from  Ref.  13.  Theoretical  contraction  ratios  for 
mitre  bends  (Ref.  14) -were  used  to  calculate  the  effect  of 
setback  ratio  on  nozzle  airflow  match  as  shown  in  Figure  37. 

The  theory  predicts  a  lower  airflow  match  than  indicated  by  test 
data,  probably  due  to  spillage  around  the  sides  of  the  plate 
for  the  test  model.  Spillage  is  not  possible  for  the  two- 
dimensional,  incompressible  potential  flow  model. 

During  Part  1C  of  this  program,  a  hinged  external  deflector  was 
tested  statically  to  determine  its  performance  both  as  a  thrust 
reverser  and  as  a  vectoring  nozzle.  Data  correlations  were 
developed  showing  the  effects  of  setback  ratio  on  airflow  match 
and  vectoring  efficiency.  The  test  results  and  data  correlations 
are  discussed  in  Volume  21. 


theoretical 
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DISCHARGE  COEFFICIENT.  CD,  OR  VELOCITY  COEFFICIENT.  Cy 


DATA:  REF  8 


0  =95° 


Figure  32:  LOBSTERTAIL  NOZZLE  PERFORMANCE  FOR 
95  DEGREE  VECTOR  ANGLE 
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DATA:  REF  8 


Figun  33:  EFFECT  OF  VECTOR  ANGLE  AND  NOZZLE 

PRESSURE  RATIO  ON  LOBSTERTAIL  NOZZLE 
PERFORMANCE 


DISCHARGE  COEFFICIENT,  Cn,  OR  VELOCITY  COEFFICIENT.  C 


DATA:  REF  12 


PLATE  LENGTH/NOZZLE  HEIGHT  RATIO 


0 


2.0 


IX) 

SETBACK  RATIO  -  X/h 


Figure  37:  MITRE  BEND  DA  TA,  CONTRACTION  COEFFICIENT 
VS  SETBA  CK  DISTANCE 
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2.1.4  Cascade  Lattice  Loss  Correlation 


> 


Correlations  were  employed  to  predict  losses  across  TR  or  TV 
cascade  lattices  in  terms  of  velocity  and  discharge  coefficients. 
The  prediction  method  employs  pressure  loss  data  of  Ainley  and 
Mathieson  (Ref.  15)  and  a  momentum  thickness  correlation 
defeloped  by  Stewart  (Ref.  16).  Ainley  and  Mathieson's  data 
are  used  to  obtain  total  pressure  loss,  exit  Mach  number  and 
exit  flow  angle.  The  entrance  and  exit  flow  properties 
used  to  calculate  reaction  across  the  blade  row.  Reaction  is 
used  in  Stewart's  correlation  for  trailing  edge  momentum 
thickness.  Velocity  and  discharge  coefficients  are  readily 
determined  from  boundary  layer  momentum  and  displacement 
thicknesses . 


A  commonly  used  system  for  defining  the  geometry  of  a  blade 
row  and  the  flow  angles  relative  to  a  blade  row  is  illustrated 
in  Figure  38.  Flow  inlet  angle  0(,  is  a  required  input  for  the 
analysis.  For  thrust  reverser  lattices,  CX,  may  vary  from  0  de¬ 
grees  to  90  degrees  along  the  blade  row.  Consequently,  a 
logical  range  of  values  should  be  used  when  analyzing  a  particu¬ 
lar  design.  Also,  note  that  the  values  of  flow  outlet  angle  0(z 
are  numerically  negative  in  Figure  38 .  For  thrust  reverser 
lattices,  discharge  angle  is  related  to  by  the  following 
equation. 


30  +■  <Xt 


(27) 


A  family  of  profile  loss  curves  are  shown  in  Figure  39  for 
reaction  blades  at  low  Mach  number  (M<-0.5),  high  Reynolds 
number  (Re  =  2  x  l(r),  and  zero  incidence.  Profile  losses  are 
presented  in  terms  of  pressure  loss  coefficient,  Yp. 


YP  = 


loss  of 

total  pressure 

r  total  pressure' 
[at  blade  outlet 

- 

Static  pressure] 
at  blade  outletj 

(28) 


A  family  of  curves  for  profile  losses  of  impulse  blades  is  shown 
in  Figure  40.  Profile  losses  of  blades  intermediate  between 
reaction  and  impulse  blades  are  interpolated  by  the  following 
equation : 


^  -WO]} 


(29) 


In  order  to  find  the  variation  of  profile  loss  with  incidence 
angle,  it  is  necessary  to  find  the  stalling  incidence  (i  ) , 
defined  as  the  incidence  at  which  the  profile  loss  is  twice  the 
minimum  loss.  The  variation  of  stalling  incidence  and  flow 
outlet  angle  with  pitch  to  chord  ratio  is  shown  in  Figure  41. 

The  effects  of  blade  inlet  angle  and  flow  outlet  angle  on  stalling 
incidence  are  shown  in  Figure  42 .  The  stalling  incidence  found 
from  Figures  41  and  42  is  used  to  find  the  relative  profile  loss 
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PROFILE  LOSS  COEFFICIENT.  Y 


I 
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Figure  40:  PROFILE  LOSSES  FOR  IMPULSE  BLADES  AT  ZERO  INCIDENCE 
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VARIATION  OF  STALLING  INCIDENCE  WITH  a  2  AND  S/C 


I 


Figure  41:  VAR  I  A  TION  OF  STALLING  INCIDENCE  AND  FLOW  OUTLET 
ANGLE  WITH  PITCH/ CHORD  RATIO 
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M  <  G.6 
Re  *  2  x  10® 


-1.0  -0.5  0  0.5  1.0 

01 

a2(S/C  *  0.75) 


Figure  42:  VARIATION  OF  STALLING  INCIDENCE  WITH  dLADE 
INLET  ANGLE  AND  FLOW  OUTLET  ANGLE 
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and  change  in  flow  outlet  angle  from  Figure  43.  For  positive 
incidence  angles ,  the  flow  outlet  angle  decreases  with  increasing 
losses.  Effects  of  exit  Mach  number  M,  are  shown  in  Figures 
44  and  45.  Flow  outlet  angle  is  obtained  from  Figure  44a  for 
exit  Mach  number  M2  0.5  and  Figure  44b  for  M.  =  1.0.  Linear 
interpolation  is  used  for  intermediate  Mach  numbers.  The  effect 
of  exit  Mach  number  and  blade  trailing  edge  curvature  on  re¬ 
lative  blade  profile  loss  is  given  in  Figure  45. 

Effects  of  Reynolds  number  are  shown  in  Figures  46  and  47.  Pro¬ 
file  losses  and  flow  outlet  angle  decrease  as  Reynolds  number 
increases. 


Entrance  and  exit  flow  properties  (Figures  39  to  47)  are  used  to 
calculate  reaction  across  the  blade  row: 


(30) 


where  V.,  and  V-  are  entrance  and  exit  velocity.  A  correlation 
between^reaction  and  trailing  edge  momentum  thickness  was 
established  by  Stewart  (Ref.  19)  as  shown  in  Figure  48.  The 
following  equation  corrects  for  Reynolds  number  effects. 


(31) 


where  (  is  obtained  from  Figure  48  and  baseline 

Re.  =300,000.  The  fraction  of  momentum  loss  based  on  actual  weight 
flow  is  given  by  the  following  equation  for  a  >.wo-dimensional 
cascade. 


(32) 


where  total  boundary  layer  displacement  thickness  at  the  blade 
trailing  edge  6^  is  related  to  momentum  thickness  by  the 
equation 


H  * 


it 


(33) 


Data  and  theory  for  form  factor  H  are  shown  in  Figure  49.  The 
term  considers  the  blades  to  have  infinite  length.  In  a 

cascade,  however,  there  are  end  walls  and  stiffeners  that  also 
have  boundary  layer  losses  that  are  appreciable.  Losses  on 
these  surfaces  are  assumed  equal  to  the  average  loss  occuring 
on  the  blade  surfaces,,  Therefore,  the  factor  to  correct  for 
end  walls  is  given  by  the  area  ratio 
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I  a,  ATi  GIVING  Yp  |  RELATIVE  PROFILE  LOSS.  Un/Y_  (i  =  0) 

*  PMIN  p  p 


VARIATION  OF  RELATIVE  PROFILE 
LOSS  WITH  RELATIVE  INCIDENCE 


Figure  43:  VAR! A  TION  OF  LOSS  AND  OUTLET  ANGLE  WITH  INCIDENCE 


BLADE  PITCH/TRAILING  EDGE  CURVATURE,  S/e 


b)  FLOW  OUTLET  ANGLE  WHEN  M2  =  1.0 

Figure  44:  VAR  I A  TION  OF  FLOW  OUTLET  ANGLE  WITH  MACH  NUMBER 
AND  BLADE  TRAILING  EDGE  CURVATURE 


RELATIVE  PROFILE  LOSS 
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Figure  46:  EFFECTS  OF  REYNOLDS  NUMBER  ON  PROFILE  LOSSES  FOR  CASCADE  LATTICES 
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,a2(Re  =  2x  lofy 


Figure  47:  EFFECTS  OF  REYNOLDS  NUMBER  ON  FLOW 
OUTLET  ANGLE  FOR  CASCADE  LATTICES 
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MOMENTUM  THICKNESS  PARAMETER, 


FORM  F 


Figure  49:  FORM  FACTOR  DA TA  FOR  CASCADE  LA  TTICES 
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For  a  cascade  lattice  with  no  stiffeners,  it  can  be  shown  that 


A so  ^  .  ,  cos  $ 

Aio  '  1  **  (35) 


where  $  is  stagger  angle  (see  Figure  38)  and  blade  aspect  ratio 
A*  =  blade  length/chord.  If  the  cascade  lattice  has  stiffeners, 
their  wetted  area  should  be  included  in  the  term  A  ^  of 


Equation  3-i.  The  area  ratio  A^/Ajp  is  used  to  correct  for 
losses  on  end  walls  and  stiffeners. 


=  I'M 


A3D 


Atp 

Velocity  coefficient  is  given  by 


(36) 


Cv  -  I  ^  UKkltrexpansioM  losses  (Fij.C) 

Discharge  coefficient  is  given  by 


(37) 


l“Hm 


2,3D 


(38) 


The  prediction  method  computes  Cv  and  C  for  cascade  lattices. 
To  obtain  reverser  efficiency,  vector  efficiency,  and  airflow 
match,  the  following  equations  are  used: 


„  _  Cv cos©r 

^  '  Cv* 


(39) 

(40) 


Cp  Am 
Af 


(41) 


2.2  Task  1. 2 — Construct  Computerized  Analytical  Models 

The  purpose  of  Task  1.2  was  to  develop  computer  programs  for 
predicting  TR  and  TV  nozzle  performance  and  evaluating  TR  and 
TV  influence  on  the  total  airplane/engine  system.  Three  computer 
programs  were  developed: 

1)  Jet  Trajectory  and  Spreading  Program  (TEM-356A) 

2)  Reingestion  Prediction  Program  (TEM-356B) 

3)  TR  and  TV  System  Performance  Program  (TEM-357) 
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The  numbers  in  parenthesis  are  permanent  identification 
numbers  assigned  to  the  programs  by  Boeing  Computer  Services, 

Inc.  Further  description  of  the  programs  is  provided  in  the 
following  sections. 

2.2.1  Jet  Trajectory  and  Spreading  Program 

The  purpose  of  the  Jet  Trajectory  and  Spreading  Program 
TEM-356A  is  to  predict  the  shape  of  the  exhaust  plume  emanating 
from  a  TR  or  TV  nozzle.  The  exhaust  plume  definition  is  used 
to  predict  potentially  severe  aerodynamic  interference  and  control 
problems  that  could  be  caused  by  plumes  impinging  on  or  passing 
close  to  flight  control  surfaces.  The  program  also  is  used  to 
provide  definition  of  TR  exhaust  plumes  to  the  Reingestion 
Prediction  Program  TEM-356B  as  described  in  Section  2.2.2. 

As  an  example  of  the  program's  usefulness  and  purpose,  Figures 
50  to  52  show  a  thrust  reverser  exhaust  plume  for  a  two-engine 
STOL  transport.  The  jet  plume  was  computed  by  TEM-356A  to 
determine  the  jet  trajectory  and  plume  shape  relative  to  the 
aerodynamic  control  surfaces  at  110  knot  freestream  velocity. 

Most  STOL  transports  are  expected  to  land  at  about  90  knots, 
so  110  knots  should  represent  a  more  critical  condition  for 
tail  interference  effects.  The  thrust  reverser  efflux  appears 
to  impinge  on  the  horizontal  stabilizer  in  the  side  view.  How¬ 
ever,  the  front  view  shows  that  the  lower  vortex  lobe  is  outboard 
of  the  stabilizer.  The  analysis  indicates  that  the  thrust  re¬ 
verser  plume  would  not  cause  any  severe  aerodynamic  interference 
problems  by  impinging  on  flight  control  surfaces. 

A  diagram  showing  inputs  and  outputs  of  the  Jet  Trajectory  and 
Spreading  Program  is  given  in  Figure  53.  Inputs  consist  of 
TR  or  TV  nozzle  position,  orientation,  exhaust  flow  direction, 
the  type  of  jet  and  nozzle  and  freesteam  flow  conditions.  Out¬ 
put  consists  of  the  jet  centerline  in  (x,  y,  z)  space  and  con¬ 
tours  of  the  jet  cross  sections  defining  the  boundaries  between 
the  TR  or  TV  exhaust  flow  and  freestream  flow.  The  analysis 
uses  an  empirical  equation  *or  the  jet  centerline,  empirical 
jet  spreading  coefficients,  and  theory  for  the  jet  cross 
section.  The  following  sections  describe  the  analysis  in  more 
detail. 

Jet  Trajectory  Equation 

There  are  numerous  empirical  equations  that  predict  the  trajectory 
of  a  jet  in  a  crossflow  (Ref.  17  to  24).  Table  IV  summarizes 
empirical  equations  for  the  jet  trajectory,  together  with  their 
respective  ranges  of  validity.  The  empiric  il  equations  were 
obtained  by  curve-fitting  data  from  flow  visualization  experi¬ 
ments.  Comprehensive  reviews  of  the  empirical  equations  were 
made  by  Filler  (Ref.  25)  and  Margason  (Ref.  24).  Filler  con¬ 
cluded  that  empirical  equations  of  Ivanov,  Shandorov,  and 
Margason  compared  favorably  within  their  respective  ranges  of 
validity.  Margason  also  concluded  that  Ivanov's  equation  provided 
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Figure  51:  STOL  TRANSPORT  THRUST  REVERSER 
PLUME  END  VIEW 


Figure  52:  STOL  TRANSPORT  THRUST  REVERSER 
PLUME  PLAN  VIEW 


JET  TRAJECTORY  AND  SPREADING  PROGRAM 


Figure  53:  JET  TRAJ ECTOR Y  AND  SPREADING 
PROGRAM  DIAGRAM 
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EQUATION  oco  qj/<Jco  I  m/t«»  reference 
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TABLE  IV  SUMMARY  CF  EMPIRICAL  EQUATIONS  FOR  JET  TRAJECTORY  AND  RANGE  OF  VALIDITY 


a  good  fit  to  his  data  at  most  of  the  deflection  angles 
between  30  and  150  degrees.  The  range  of  validity  for 
Ivanov's  equation  is  clearly  known  and  closely  meets 
conditions  of  thrust  reverser/vectoring,  namely 
0  £  qj/q^  —  200  and  30  degrees^  ec-g  150  degrees.  For 

these  reasons,  Ivanov's  equation  was  selected  for  use  in  the 
Jet  Trajectory  and  Spreading  Program.  Trajectories  comput¬ 
ed  using  Ivanov's  equation  are  compared  to  test  data  in 
Figures  54  to  56. 

Ivanov  also  experimented  with  rectangular  jets  (1:5  £=a/b£:  5:1) 
and  concluded  that  to  a  first  approximation  the  equation  for 
a  circular  jet  is  adequate  provided  the  hydraulic  diameter 
is  used  instead  of  the  initial  jet  diameter. 

For  two-dimensional  jets,  the  equation  of  Vizel  and  Mostinskii 
(Ref.  23)  was  selected. 
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Experimental  data  for  jet  penetration  coefficient  C  are  shown 
in  Figure  57  from  Ref.  26.  The  straight  lines  drawn  through 
the  data  are  used  in  program  TEM-356A.  The  line  for  <X  =  90 
degrees  is  used  for  discharge  angles  less  than  90  degrees 
and  the  line  for  o (  =  135  degrees  is  used  for  angles  greater 

than  135  degrees.  Interpolation  is  used  to  obtain  jet 
penetration  coefficients  for  90  degrees^  or  -c  135  degrees. 
Comparisons  of  Vizel  and  Mostinskii 's  results  with  test  data 
are  shown  in  Figures  58  and  59. 

Jet  Spreading  Coefficients 

Experimental  jet  spreading  coefficients  obtained  for  a 
round  jet  discharging  at  90  degrees  to  a  freestream  flow 
(Ref.  22)  are  used  to  predict  thickness  and  width  of  the 
exhaust  plume.  A  lack  of  data  exists  for  other  discharge 
angles  and  jet  shapes.  Moreover,  there  is  considerable 
disagreement  for  jet  thickness  to  width  ratio,  S/w,  between 
the  data  of  several  investigators. 

Pratte  and  Baines'  jet  thickness  spreading  coefficient  data 
for  a  round  jet  in  crossflow  are  presented  in  Figure  60  from 
Ref.  22.  The  following  empirical  equations  provide  a  good 
curve  fit  to  the  data. 
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ttgun  54:  COMPARISON  OF  JET  TRAJET&.  <  EQUATION  WITH  DA  TA  FOR  an  =  135  DEGREES 
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JET  PENETRATION  COEFFICIENT.  Cx 


VELOCITY  RATIO,  Vj/Uoo 


Figure  57:  JET  PENETRATION  COEFFICIENT  DATA  FOR 

VIZEL  AND  MOSTiNSKII'S  TRAJECTORY  EQUATION 
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These  empirical  equations  are  used  downstream  of  the  jet's 
potential  core.  In  the  potential  core  region,  the  thickness 
to  diameter  ratio  is  equal  to  1.0.  Setting  S/d  “  1.0  in 
Equation  43  and  solving  for  s/d  gives 


'  d©/  tote* full 


O.IS93 
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core 

Summarizing,  the  following  equations  are  used  to  compute  jet 
thickness . 
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Pratte  and  Baines  recommend  a  jet  width-to-thickness  ratio  of 
1.4  based  upon  their  flow  visualization  experiments.  This 
value  compares  to  width/thickness  ratio  w/S  st  1  measured 
by  Gerend  (Ref.  28),  w/f  ■  4  recommended  by  Wooler  (Ref.  29) 
and  w/S  *  5  suggested  by  Abramovich  (Ref.  17).  Pratte  and 
Baines'  w/S  ■  1.4  was  selected  because  the  test  section  was 
considered  large  enough  to  represent  full-scale  jet  spread¬ 
ing  characteristics  and  because  of  superior  data  repeatability. 

In  the  potential  core  region,  jet  width  is  obtained  by  inter¬ 
polating  between  w/S  «  1  at  the  origin  and  w/S  ■  1.4  at 
the  end  of  the  potential  core.  Equations  for  jet  width  are 
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Jet  Cross-Section  Analysis 

A  number  of  recent  theoretical  investigations  have  been  made 
using  a  theory  developed  by  Chang  (Ref.  30  to  33).  Chang 
developed  a  two-dimensional  potential  flow  theory  that  predicts 
the  vortex  rollup  of  a  cylinder  in  crossflow.  Results  of 
Chang’s  analysis  are  shown  in  Figure  61.  The  cross  section  is 
defined  by  48  points  and  was  computed  for  a  velocity  ratio 
vj/uoo  *  8*  Rollup  of  the  cylinder  into  two  counter¬ 
rotating  vortices  is  evident.  It  should  be  noted  that  Chang's 
theory  predicts  the  cross-sectional  shape  of  the  jet,  but 
not  the  trajectory.  The  cross  sections  are  placed  on  the 
trajectory  predicted  by  Ivanov's  equation. 


Cross  sections  predicted  by  Chang's  theory  are  remarkably 
similar  to  vortex  rollup  phenomena  observed  in  flow 
visualization  experiments.  Photographs  showing  vortex 
rollup  are  shown  in  Figure  62.  A  brief  discussion  of 
Chang's  theory  is  given  in  Appendix  I. 


Jet  Plume  Construction  Procedure 


The  following  general  procedure  is  used  to  construct  the 

jet  plume.  The  frees tream  velocity  vector  (defined  by 

angles  <X  and  ?  )  and  the  jet  initial  direction  (defined  by 

vector  components  t  ,  t  ,  t  )  define  a  plane,  as  shown  in 

**  y  z 

the  following  sketch. 
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ISOMETRIC  VIEW 


The  jet  centerline  is  assumed  to  lie  in  the  plane  defined  by 
vectors  and  T. 

Jet  cross  sections  are  computed  for  each  station  along  the 
jet  centerline  using  Chang's  analysis  for  round  or  rectangular 
jets.  Thickness  and  width  spreading  coefficients  are  used 
to  define  gross  deformation  of  the  jet  as  shown  in  the  following 
sketch. 


The  program  then  centers  the  cross  section  on  the  axis  and 
scales  the  jet  cross-section  to  the  correct  size,  as  shown 
in  the  following  sketch. 
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Geometric  Capability 

The  Jet  Trajectory  and  Spreading  Program  was  formulated 
with  considerable  flexibility  so  that  it  can  be  applied  to 
a  wide  vr.riety  of  TR/TV  geometries.  The  program  predicts 
four  types  of  TR/TV  exhaust  plumes.  Inputs  and  calculation 
procedures  unique  to  each  jet  type  are  discussed  below. 

1)  TYPE  =  1.  Circular  jet  cross-section. 

This  type  of  jet  models  the  flow  from  round  nozzles, 
e.g.,  three-bearing  nozzle,  spherical  eyeball,  or  a  pure 
lift  engine  buried  in  a  fuselage  or  pod.  Inputs  for 
TYPE  =  1.  are  depicted  graphically  in  Figure  63  and 
summarized  in  Table  V.  The  calculations  use  Ivanov's 
trajectory  equation  and  Pratte  and  Baines'  spreading 
coefficients.  Vortex  strengths  for  Chang's  analysis 
are  computed  using  Equation  11  in  Appendix  I.  A 
sample  case  is  given  in  Appendix  II. 

2)  TYPE  =  2.  Rectangular  jet  cross-section. 

This  type  of  jet  is  best  suited  for  cascade  or  target 
thrust  reverser  plumes  that  can  be  approximated  by  a 
rectangular  initial  cross-section.  Inputs  for  TYPE  =  2. 
are  depicted  in  Figure  64  and  summarized  in  Table  VI. 

The  calculations  use  a  hydraulic  diameter  in  Ivanov's 
trajectory  equation 

,  _  4(  »r««)  _  Z  a  b 

p« A4  b  (55) 

where  a  and  b  are  the  length  and  width  of  the  initial 
cross-section.  Pratte  and  Baines'  empirical  jet-spreading 
coefficients  are  used  for  jet  thickness  and  width. 

Vortex  strengths  for  Chang's  analysis  are  computed  from 
table  lookup  as  described  in  Appendix  I.  A  sample  case 
is  given  in  Appendix  II. 

3)  TYPE  *  2.  Two-dimensional  jet  cross-section. 

This  type  of  jet  models  the  flow  from  a  slot  nozzle  or 
high  aspect  ratio  rectangular  nozzle.  Inputs  for 
TYPE  •  3.  are  depicted  in  Figure  65  and  summarized  in 
Table  VII.  Vizel  and  Kostinskii's  quadratic  equation  is 
used  for  the  jet  centerline.  Pratte  and  Baines'  spreading 
is  used  for  the  jet  width.  Chang's  analysis  for  the 
jet  cross-section  is  not  used  because  it  was  developed 
for  a  discrete  jet  shape,  not  an  infinite  jet  sheet.  A 
sample  case  is  given  in  Appendix  II . 

4)  TYPE  =  4.  Annular  jet  cross-section. 
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TABLE  VI  INPUTS  FOR  TYPE  =  2.,  RECTANGULAR  JET  CROSS  SECTION 


length.  Any  consistent  length  scale  is  acceptable,  e.g. ,  inches,  fee 
ND  indicates  non-dimensional. 


ISOMETRIC  VIEW 


TABLE  VII  INPUTS  FOR  TYPE  =  3.,  TWO-DIMENSIONAL  JET  CROSS  SECTION 
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NOTE:  *L  indicates  length.  Any  consistent  length  scale  is  acceptable,  e.g.,  inches,  feet 

or  meters.  ND  indicates  non-dimensional. 


This  type  of  jet  models  the  flow  from  annular  target 
or  cascade  thrust  reversers.  Inputs  for  TYPE  *  4. 
are  depicted  in  Figure  66  and  summarized  in  Table  VIII. 

The  calculations  use  Vizel  and  Mostinskii's  trajectory 
equation  and  Pratte  and  Baines'  thickness  spreading 
coefficient.  Chang's  analysis  is  not  applicable  for 
annular  shaped  jets. 

As  shown  in  Figure  66,  the  output  consists  of  the  plume 
front  surface,  centerline,  and  rear  surface  at  a  number 
or  radial  cutting  planes.  Origins  of  the  jet  arms  lie 
in  a  plane  perpendicular  to  the  &  vector.  Their 
position  in  that  plane  is  defined  by  initial  radius  R 
and  clock  angle  6. 

The  four  types  of  jets  offer  sufficient  geometric  capability 
for  most  thrust  reverser  and  vectoring  nozzle  designs. 

Each  type  of  jet  requires  certain  inputs  describing  the  jet's 
position  and  orientation  in  space.  All  inputs  and  outputs 
are  made  in  an  airplane  (x,y,z,)  reference  coordinate  system 
defined  in  Figure  63.  The  x  axis  is  parallel  to  the  body  of 
the  airplane,  the  y  axis  is  in  the  spanwise  direction  of  the 
right  half  of  the  wing,  and  the  z  axis  is  upward.  The  origin 
is  chosen  at  any  convenient  point,  usually  at  the  airplane 
nose.  The  angles  of  attack  and  yaw  are  defined  in  terms  of 
the  direction  of  the  freestream  with  respect  to  the  reference 
coordinate  system,  as  shown  by  the  sketch  on  page  85. 

2.2.2  Reingestion  Prediction  Program 

Reingestion  occurs  when  the  deflected  TR  and  TV  exhaust  flow 
penetrates  forward  into  the  freestream  flow  and  is  captured 
by  the  inlet  flow.  Pressure  and  temperature  distortions  caused 
by  reingestion  of  the  exhaust  flow  reduce  engine  thrust  and 
compressor  stall  margin  and  can  lead  to  engine  surge. 

The  purpose  of  the  Reingestion  Prediction  Program  is  to 
predict  the  occurrence  of  reingestion  for  arbitrary  thrust 
reverser  and  airplane  configurations  as  a  function  of  geometry 
and  flow  conditions.  The  program  is  intended  for  use  as  a 
diagnostic  tool  to  predict  severe  reingestion  problems.  The 
program  does  not  predict  inlet  temperature  rise  or  the  relative 
severity  of  reingestion.  However,  the  program  printout  is 
designed  such  that  in  many  cases  judgment  can  be  made  about 
the  relative  severity  of  reingestion. 

In  order  to  develop  reingestion  prediction  methods,  reingestion 
was  classified  into  four  categories: 
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TABLE  VIII  INPUTS  FOR  TYPE  =  4.,  ANNULAR  JET  CROSS  SECTION 
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NOTE:  *  L  indicates  length.  Any  consistent  length  scale  is  acceptable,  e.g.,  inches,  feet 

or  meters.  ND  indicates  non-dimensional. 


o  Crossflow  reingestion 
o  Self  or  closed  loop  reingestion 
o  Near-field  fountain  reingestion 
o  Far-field  fountain  reingestion 

The  computer  program  treats  each  category  of  reingestion 
separately.  Computer  program  inputs  are  the  TR/TV  and  inlet 
geometry  and  flow  conditions  as  diagrammed  in  Figure  67. 
Output  consists  of  diagnostic  statements  predicting  whether 
or  not  reingestion  occurs  for  the  case  in  question.  Methods 
used  to  predict  each  category  of  reingestion  are  discussed 
below. 

Crossflow  Reingestion 

Crossflow  reingestion  occurs  when  the  exhaust  flow  from  one 
engine  is  captured  by  the  inlet  flow  of  an  adjacent  engine. 
The  most  common  occurrence  of  crossflow  reingestion  is  in 
the  outboard  engine  of  a  four-engine  airplane,  as  shown  in 
Figure  68. 

The  method  used  to  predict  crossflow  reingestion  employs 
thrust  reverser  exhaust  plumes  calculated  by  the  Jet 
Trajectory  and  Spreading  Program.  The  exhaust  plumes  are 
checked  for  intersections  with  inlet  streamtubes.  The 
program  prints  out  the  number  of  points  intersecting  an 
inlet  streamtube  at  each  jet  station. 

To  predict  crossflow  reingestion,  definition  of  the  inlet 
streamtube  is  required.  Streamtubes  were  calculated  for  a 
representative  STOL  transport  inlet  using  axisymmetric  com¬ 
pressible  potential  flow  program  TEM-095  (Ref.  34).  The 
inlet  hilite  to  throat  contraction  ratio  is  approximately 
1.3,  diffuser  angle  equals  5  degrees,  and  throat  Mach  is  0.6. 
Computer-generated  plots  of  streamlines  and  symbols  indicat¬ 
ing  Mach  contours  are  shown  in  Figures  69  and  70  for  free- 
stream  velocities  of  2  and  10  knots,  respectively.  For  low 
freestream  velocities,  the  stagnation  streamline  attaches 
on  the  cowl  outer  surface  far  downstream  of  the  cowl  leading 
edge. 

The  stagnation  streamline  defines  the  pre-entry  streamtube 
used  in  predicting  crossflow  reingestion.  Streamtubes  are 
summarized  in  Figure  71  for  inlet  hilite  to  freestream 
velocity  ratio  from  137.3  to  2.5.  Streamtubes  for  other 
velocity  ratios  are  obtained  by  quadratic  interpolation  or 
extrapolation  to  obtain  the  streamtube  radius  R  at  any 
station. 

Arbitrary  crosswind  angle  is  accounted  for  by  aligning  the 
axisymmetric  streamtube  in  the  direction  of  the  freestream 
flow,  as  shown  in  Figure  72.  A  local  inlet  coordinate  system 
is  constructed  along  the  streamtube  axis.  The  (x,y,z)  co¬ 
ordinates  of  the  exhaust  plume  computed  by  the  Jet  Trajectory 
and  Spreading  Program  are  transferred  into  this  inlet  co¬ 
ordinate  system  and  checked  for  possible  intersections. 
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Figure  68:  CROSS  FLOW  REINGESTION 
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Figure  69:  STREAMLINES  AND  MACH  CONTOURS  FOR  STOL 
TRANSPORT  INLET 
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Figure  70:  STREAM L INES  AND  MACH  CONTOURS  FOR  STOL 
TRANSPORT  INLET 
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Existing  jet  penetration  data  (Ref.  35)  show  that  the  thrust 
reverser  plume  is  extremely  unsteady  and  that  large  turbulent 
fluctuations  occur  in  the  mixing  region  with  the  freestream 
flow.  To  account  for  flowfield  unsteadiness  in  predicting 
crossflow  reingestion,  the  correlation  shown  in  Figure  73 
is  used.  Data  are  shown  for  the  jet  maximum  penetration  point 
(MPP)  as  a  function  of  velocity  ratio  and  discharge  angle. 

The  scatter  radius  is  given  by 


/-v  /  va**r  .©.fctsi 


From  the  nozzle  exit  to  the  MPP,  a  linear  interpolation  is 
used  to  obtain  scatter  radius. 


rs  _  s/d  (  r>\ 
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Downstream  of  the  maximum  penetration  point,  the  scatter 
radius  at  the  MPP  is  used 


The  scatter  radius  is  a  measure  of  uncertainty  in  predicting 
the  plume's  position  due  to  large-scale  turbulent  fluctuations 
of  exhaust  flow.  Because  of  this  uncertainty,  the  plume's 
forward  boundary  is  increased  by  the  scatter  radius.  The 
adjusted  plume  surface  is  checked  for  crossflow  reingestion. 

Crossflow  reingestion  results  are  shown  in  Figure  74.  The 
configuration  is  an  early  version  of  the  747  airplane  with  long 
duct  nacelles  and  annular  target  thrust  reversers.  The  con¬ 
figuration  was  tested  for  reingestion  characteristics  in  a 
low-speed  wind  tunnel  test.  The  computer  program  predicts 
crossflow  reingestion  at  about  90  knots.  The  wind  tunnel 
test  indicated  crossflow  reingestion  at  90  knots  (Ref.  36). 

The  main  limitation  to  the  crossflow  prediction  method  is 
that  inlet  and  exhaust  flows  are  assumed  to  act  independently. 
Effects  of  inlet  suction  on  the  exhaust  flow  are  not  included. 
Mutual  interactions  between  impinging  jet  flows  are  also 
not  included. 

Self  or  Closed  Loop  Reingestion 


Self  or  closed  loop  reingestion  occurs  when  the  TR  or  TV 
exhaust  flow  forms  a  continous  closed  path  back  into  the 
engine  inlet.  Two  types  of  self-reingestion  are  shown  in 
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an 


0  10  20  30  40  SO  60 

V:  0.94  A  CQC 

2.97  (— 1— )  (1  -0.734  SIN  Ob8&0) 

Figun  73:  JET  PENETRA  TION  COE  RE  LA  TION 

ILLUSTRATING  DATA  SCATTER  DUE  TO 
TURBULENT  FLUCTUA TIONS  OF  FLOW 
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CROSSFLOW  REINGESTION 


Figurt  74:  CROSSFLOW  REINGESTION  RESUL  TS  FOR 
90  KNOTS  FREESTREAM  VELOCITY 


Figure  75.  The  first  type  is  caused  by  the  thrust  reverser 
exhaust  flow  attaching  to  the  cowl  surface  and  penetrating 
forward  into  the  oncoming  flow,  where  it  is  captured  by  the 
inlet  flow.  For  the  second  type,  the  exhaust  flow  does  not 
attach  to  the  cowl,  but  the  proximity  of  inlet  and  nozzle 
causes  exhaust  flow  to  be  captured.  This  type  of  self- 
reingestion  is  treated  by  the  procedure  discussed  previously 
for  crossflow  reingestion. 

In  order  to  predict  self-reingestion  due  to  exhaust  flow  re¬ 
attachment,  it  was  first  necessary  to  develop  reattachment 
criteria.  For  a  two-dimensional  jet  discharging  obliquely 
from  a  plate,  it  is  well  known  that  under  certain  conditions 
the  jet  bends  towards  the  plate  and  reattaches  due  to  re¬ 
duced  pressures  caused  by  entrainment  in  the  separation 
bubble.  As  shown  by  the  data  in  Figure  76,  the  jet  will 
always  be  separated  for  discharge  angles  greater  than  about  60 
degrees.  As  the  discharge  angle  decreases,  the  value  of 
plate  length  over  slot  width  required  for  attached  flow  also 
decreases.  Figure  76  is  used  to  predict  reattachment  for 
two-dimensional  and  annular  jets  (TYPE  *=  3.  and  4.). 

The  reattachment  of  a  three-dimensional  and  annular  jet  behaves 
differently  from  a  two-dimensional  jet.  Freestream  flow 
enters  the  separation  bubble  from  the  side  as  illustrated  in 
Figure  77,  thus  relieving  reduced  pressures  in  the  bubble. 
Presence  of  a  counterflow  further  complicates  the  flow  field. 
Considerable  effort  was  made  in  searching  for  data  or  theory 
to  predict  the  reattachment  of  a  three-dimensional  jet.  No 
theory  and  very  little  data  were  discovered.  Test  data 
from  an  SST  thrust  reverser  concept  with  a  30-degree  dis¬ 
charge  angle  revealed  a  strong  reattachment  under  static 
conditions.  A  schematic  of  the  thrust  reverser  and  flow 
field  is  shown  in  Figure  78a.  Oil  flow  visualization 
proved  that  reattachment  occurred  along  the  centerline  of  the 
second  door  as  shown  in  Figure  78b. 

Reingestion  testing  conducted  during  the  development  of  the 
707  cascade  thrust  reverser  resulted  in  the  following  general 
criteria  for  jet  discharge  angles 

f  55  degrees,  jet  will  no  reattach 
0  *  ■<  45  degrees,  jet  may  or  may  not  reattach  (60) 

t  30  degrees,  jet  will  reattach 

Based  on  the  707  and  SST  test  data,  the  criteria  given  by 
Equation  (60)  are  used  to  predict  reattachment  for  discrete 
jets  (TYPE  *  1.  and  2.). 

No  method  has  been  developed  to  predict  whether  or  not  re¬ 
ingestion  occurs  once  a  jet  has  reattached  to  the  cowl.  Con¬ 
ceivably,  it  could  separate  again  in  presence  of  a  strong 
counter flow.  Experience  has  shown,  however,  that  reingestion 
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a)  SELF  REINGESTION  DUE  TO  THRUST  REVERSER  EXHAUST  FLOW 
ATTACHMENT  ON  COWL  SURFACE 


b)  SELF  REINGESTION  DUE  TO  PROXIMITY  OF  INLET  AND  THRUST 
REVERSER  NOZZLE 


Figun  75:  SELF  REINGESTION 


PLATE  LENGTH 
SLOT  WIDTH 


SYMBOL  REFERENCE 


0  BOURQUE  AND  NEWMAN,  REF  37 
0  NACA  TN  4272,  NACA  TN  4377,  REF  38 


DISCHARGE  ANGLE,  0  ,  DEGREES 


Figure  76:  EFFECT  OF  DISCHARGE  ANGLE  AND  PLA  TE  LENGTH 
ON  THE  REA  TTACHMENT  OF  A  TWO  DIMENSIONAL 
INCOMPRESSIBLE  JET 
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SEPARA' 

BUBBLE 


b)  OIL  FLOW  VISUALIZATION 


Figure  78: 


EXHAUST  FLOW  REA  TTACHMENT  ON  SST  THRUST 
REVERSER 


usually  occurs  if  the  flow  reattaches.  Consequently,  the 
approach  taken  for  self-reingestion  is  that  the  reattachment 
criteria  also  predict  reingestion.  Thus,  if  reattachment  is 
predicted  from  Figure  76  or  Equation  60,  then  reingestion 
also  is  predicted. 

There  are  several  limitations  to  the  self-reingestion  pre¬ 
diction  method.  Inlet  and  exhaust  flows  are  assumed  to  act 
independently.  Effects  of  inlet  suction  on  exhaust  flow  are 
not  included.  Secondly,  there  is  no  definitive  reattachment 
criterion  for  discrete  jets.  Also,  there  is  no  method  to  pre¬ 
dict  whether  a  reattached  jet  will  continue  to  pro¬ 
pagate  forward  or  separate  and  be  blown  rearward.  Neverthe¬ 
less,  the  method  fulfills  the  major  objective  of  providing  a 
diagnostic  tool  for  reingestion  problems. 

Near-Field  Fountain  Reingestion 

Near-field  fountain  reingestion  occurs  when  two  or  more  ad¬ 
jacent  jets  impinge  nearly  vertically  on  the  ground  plane. 
Spreading  ground  flows  collide  and  cause  "fountains"  of  high 
velocity  air  to  rise  vertically  as  shown  in  Figure  78.  These 
fountains  can  be  ingested  into  nearby  inlets.  The  prediction 
method  for  near-field  fountain  reingestion  employs  existing 
NASA  test  data  from  Ref.  39.  The  test  evaluated  recircul¬ 
ation  effects  resulting  from  a  pair  of  heated  jets  imping¬ 
ing  on  a  ground  plane.  Significant  conclusions  reached 
about  conditions  causing  near-field  fountain  reingestion 
are  summarized  below: 

1)  Upwash  and  inlet  temperature  rises  were  extremely 

sensitive  to  small  nozzle  cant  angles  that  result  in 
configuration  asymmetries  as  viewed  looking  normal  to 
the  common  plane  of  the  nozzles.  A  sketch  of  the  flow- 
field  ’.,ith  canted  nozzles  is  shown  in  Figure  80a.  Con¬ 
versely,  the  upwash  and  inlet  temperature  rises  were 
relatively  insensitive  to  small  nozzle  cant  angles  that 
maintained  symmetry  as  viewed  looking  normal  to  the 
common  plane  of  the  nozzles  as  shown  in  Figure  80b. 

2}  The  inlet  flow  rate  had  almost  no  effect  on  inlet  temper¬ 
ature  rise.  The  upwash  flow  was  established  predomin¬ 
antly  by  the  jets  with  the  inlets  merely  swallowing  air 
in  their  proximity.  The  upwash  between  nozzles  generally 
was  concentrated  in  a  region  about  three  diameters  in 
width,  with  velocities  frequently  in  excess  of  .100  feet 
per  second.  No  data  for  upwash  temperature  or  velocity 
were  presented  for  the  plane  of  symmetry  between  the 
two  nozzles.  However,  flow  visualization  photographs 
and  pressure  surveys  on  the  ground  plane  indicate  the 
fountain  was  at  least  six  diameters  wide  along  the  stag¬ 
nation  streamline  dividing  flows.  It  can  therefore  be 
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79:  NEAR-FIELD 


b)  NOMINAL  FLOW  FIELD  WITH  SYMMETRIC  NOZZLES  AND 
EQUAL  NOZZLE  PRESSURE 

Figure  80:  NEAR-FIELD  FOUNTAIN  FLOW  FIELDS 


assumed  that  inlets  located  parallel  to  the  symmetry 
plane  also  reingest  flow  from  the  fountain. 

■t* 

3)  The  jets  merged  if  they  were  located  close  together  and 
high  off  the  ground.  The  criterion  for  a  pair  of  parallel 
jets  to  merge  is  given  by 

(Dj.*  <61>  : 

Wlftrje 

where  <x  is  the  half  angle  of  the  spreading  jet  plume ,  S 
is  the  spacing  between  jets,  and  H  is  the  height.  Using 
a  spreading  angle  of  7  degrees  gives 

When  the  jets  merged,  no  fountain  or  reingestion  occurred. 

3 

4)  The  effect  of  wind  in  altering  the  near  flowfield  was 
negligible  for  wind  speeds  less  than  6  feet  per  second. 

With  the  wind  parallel  to  the  common  plane  of  the  nozzles, 
further  increases  in  wind  speed  shifted  the  fountain 
downstream,  until  at  about  16  feet  per  second,  the  up- 

wash  was  shifted  beyond  the  downstream  jet,  thus  elimin-  'j 
ating  the  fountain  between  the  jets. 

With  the  wind  at  90  degrees  to  the  common  plane  of  the 

nozzles,  no  significant  change  in  either  the  upwash 

flow  or  inlet  temperature  rise  was  detected  within  the 

range  of  wind  speeds  investigated,  i.e.,  wind  speed  <16  J 

feet  per  second.  However,  analysis  of  the  jet  trajectory 

of  the  fountain  shows  that  for  wind  speeds  greater  than 

about  25  feet  per  second,  the  fountain  will  be  swept 

downstream. 

Based  upon  these  conclusions,  certain  criteria  have  been  estab¬ 
lished  for  near-field  fountain  reinqestion  to  occur,  namely: 

1)  Near-field  fountain  reingestion  will  not  occur  for  free- 
stream  velocities  greater  than  about  25  feet  per  second 
because  the  fountain  is  swept  downstream  by  the  freestream 

flow.  I 

2)  There  must  be  a  minimum  of  two  jets  directed  downward 
for  a  fountain  to  form. 

3)  The  jet  pairs  are  checked  to  see  whether  they  merge.  If 

the  jets  merge  before  impinging  on  the  ground,  no  fountain  ' 

or  reingestion  occurs. 
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4)  If  the  cant  angles  are  symmetric  when  viewed  in  their 

common  plane#  near-field  fountain  reingestion  will  occur 
for  cant  angles  up  to  ±  20  degrees.  If  the  cant  angles 
are  not  symmetric#  near-field  reingestion  will  occur 
for  cant  angles  up  to  about  6  degrees. 

The  Reingestion  Prediction  Program  checks  all  jets  against 
the  above  criteria.  If  a  pair  of  jets  satisfying  the  criteria 
are  found,  then  reingestion  diagnostic  comments  are  printed 
out. 


The  principal  limitation  of  the  near-field  fountain  reingestion 
prediction  method  is  the  fact  that  the  test  hardware  represented 
a  lift  engine  having  an  inlet  located  very  close  to  the  exhaust 
nozzle  (Figure  80).  STOL  transports  employing  thrust  vector¬ 
ing  of  the  cruise  engines  will  have  inlets  located  several 
nozzle  diameters  in  frontof  the  nozzle.  This  type  of  con¬ 
figuration  will  be  less  susceptible  to  near-field  reingestion 
than  the  NASA  hardware.  Therefore#  reingestion  predictions 
using  NASA  data  will  be  conservative. 

Far-Field  Fountain  Reingestion 

Far-field  fountain  reingestion  occurs  when  the  spreading  ground 
flow  from  one  of  more  impinging  jets  separates  from  <he  ground# 
rises  in  a  cloud,  and  is  blown  back  into  the  inlet.  Separa¬ 
tion  is  caused  by  buoyancy  forces  in  the  absence  of  a  free- 
stream  flow  or  by  interaction  with  a  headwind. 

Figure  81  depicts  the  major  elements  used  to  predict  far- 
field  reingestion.  The  dividing  streamline  between  jet  and 
freestream  flow  is  determined  first.  Secondly,  the  upper 
surface  of  the  exhaust  cloud  is  located.  Finally#  the  exhaust 
cloud  is  checked  for  intersections  with  the  inlet  streamtubes.  De¬ 
tails  of  this  procedure  are  discussed  in  the  following  para¬ 
graphs  . 

To  predict  far-field  reingestion,  the  separation  point  of  the 
spreading  ground  jet  is  required.  Separation  data  from  Ref. 

40  to  43  are  compared  in  Figure  82.  Abbott's  data  agree  fairly 
well  with  Binion's  data  and  with  NASA's  single  and  dual  pod 
data  over  a  limited  range  of  dynamic  pressure  ratios.  Dis¬ 
crepancies  'are  probably) due  to  configuration  differences  and 
measurement  inaccuracies.  Abbott  correlated  data  for  a  wide 
range  of  dynamic  pressure  ratios  and  impingement  angles  and 
developed  a  criterion  to  predict  the  separation  distance 
(Ref.  40).  Abbott *8  separation  criterion#  depicted  in 
Figure  83#  states  that  separation  occurs  at  a  distance 
where  the ’dynamic  pressure  for  a  stationary  nozzle  is 
equal  to  four  times  the  dynamic  pressure  of  the  still  air 
relative  to  a  moving  nozzle#  i.e.# 


*»%«(• 


(63) 
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Figure  81:  FAR-FiELD  FOUNTAIN  REINGESTION 
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Figure  83:  ABBOTT'S  CRITERION  FOR  PREDICTING 
FIELD  FOUNTAIN  FLOW  SEPARATION 


Abbott  conducted  extensive  dynamic  pressure  surveys  for 
stationary  nozzles  and  used  flow  visualization  for  moving 
nozzles  in  order  to  develop  this  criterion.  Results  are 
shown  in  Figure  84  for  impingement  angles  from  40  to  150 
degrees  and  dynamic  pressure  ratios  from  4  to  40,000. 

Abbott's  correlation  for  separation  distance  of  t’ie  spreading 
ground  jet  applies  only  in  the  vertical  plane  of  symmetry. 

The  general  case  of  far-field  reingestion  involves  inlets 
located  to  the  right  or  left  of  this  plane.  Therefore,  it 
is  necessary  to  predict  the  position  of  the  dividing  stream¬ 
line  out  of  the  plane  of  symmetry.  Potential  flow  theory  was 
used  by  Colin  with  good  results  to  predict  the  peeling  line 
for  a  jet  impinging  at  90  degrees  (Ref.  41).  Comparison  of 
Colin's  theory  to  test  data  is  shown  in  Figure  85.  His 
potential  flow  model  consisted  of  a  line  source  combined  with 
a  freestream  flow.  The  source  strength  was  sized  to  make 
the  dividing  streamline  match  data  in  the  plane  of  symmetry. 

Colin's  theory  was  extended  to  jets  impinging  at  arbitrary 
angles  by  specifying  a  nonuniform  outflow  velocity  distri¬ 
bution  around  the  periphery  of  the  jet's  cross  ection  pro¬ 
jected  on  the  ground  plane.  The  outflow  velocity  was 
determined  by  balancing  the  momentum  of  radial  streamlines 
leaving  the  stagnation  point.  Computer-plotted  results  are 
shown  in  Figures  86  to  88  for  0,  45,  and  90-degree  wind 
direction.  Streamlines  are  shown  emanating  from  the  jet 
stagnation  point  and  in  the  freestream  surrounding  the  zone  of 
jet  flow.  Streamlines  without  symbols  were  traced  from 
separate  computer  submissions.  Jet  impingement  angles  of 
30,  60,  90,  120,  and  150  degrees  were  run  with  0,  45,  and 
90-degree  wind  direction  at  velocity  ratios  Vj  /U*,  =  1  to  100. 
Wind  direction  has  a  pronounced  influence  on  the  dividing 
streamline  as  shown  by  Figure  89.  However,  when  the 
dividing  streamlines  are  aligned  along  the  direction  of 
approach  flow  as  in  Figure  90  the  effect  of  wind  direction 
is  seen  to  be  small.  As  expected,  velocity  ratio  had  a 
large  influence  on  the  dividing  streamline  as  shown  in 
Figure  91.  However,  the  influence  of  jet  impingement 
angle  on  dividing  streamline  shape  was  unexpectedly  small, 
indicating  that  total  jet  strength  is  more  important  in 
determining  dividing  streamline  shape  than  the  relative 
distribution  of  outflow  velocity. 

As  discussed  previously,  Abbott's  criterion  is  used  to  predict 
the  separation  point  in  the  plane  of  symmetry.  The  peeling 
line  radius  out  of  the  plane  of  symmetry  is  determined  from 
the  dividing  streamline  for  90  degree  impingement.  Because 
the  jet  impingement  angle  has  only  a  small  influence  on 
dividing  streamline  shape,  the  dividing  streamline  for  90 
degree  impingement  is  used  for  all  impingement  angles. 
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Figure  85:  DIVIDING  STREAMLINE  COMPARISON 
BETWEEN  THEORY  AND  EXPERIMENT 
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LATERAL  DISTANCE  FROM  PROJECTED  IMPINGEMENT  POINT  <Y) 


-2.000 L. 
2.00 


\  SEPARATION  POINT  \  STAGNATION  POINT 
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Figure  86:  GROUND  PLANE  STREAMLINES  FOR  3CP 
JET  IMPINGEMENT  ANGLE  AND  VELOCITY 
RATIO  V./U„=  1.0 
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DISTANCE  FROM  PROJECTED  IMPINGEMENT  POINT  (X) 

Figure  88:  GROUND  PLANE  STREAMLINES  FOR  3CP 
JET  IMPINGEMENT  ANGLE  AND  9CP  WIND 
DIRECTION 
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Figure  89:  EFFECT  OF  WIND  DIRECTION  ON  Dl  VIDING 
STREAMLINE  SHAPE 
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Figun  90:  EFFECT  OF  WIND  DIRECTION  ON  DIVIDING 
STREAMLINE  SHAPE 
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Once  the  position  of  the  dividing  streamline  is  determined, 
it  is  necessary  to  predict  the  height  of  the  exhaust  cloud. 
Abbott's  data  shown  in  Figure  92 'are  used  for  this  purpose. 

Data  scatter  is  due  in  part  to  turbulent  fluctuations 
present  in  the  exhaust  cloud.  The  top  line  is  used  to 
provide  a  factor  of  conservatism  in  the  prediction  method. 

The  radial  position  of  the  top  of  the  exhaust  cloud  is 
assumed  to  lie  above  the  vortex  center  given  by  Colin's 
data,  as  shown  in  Figure  93. 

The  final  step  necessary  to  predict  far-field  reingestion 
is  to  check  the  exhaust  cloud  upper  surface  for  possible 
intersections  with  inlet  streamtube  surfaces.  The  Re¬ 
ingestion  Prediction  Program  checks  a  limited  number  of  points 
in  the  cloud  for  intersections  with  inlet  streamtube  surfaces. 
As  shown  in  Figure  94,  points  are  checked  at  azimuth  angles 
from  30  to  330  degrees.  This  simplified  intersection  pro¬ 
cedure  was  chosen  because  the  exhaust  cloud  surface  is  not 
known  with  sufficient  accuracy  to  warrant  a  more  sophisticated 
surface  intersection  technique. 

Steps  used  to  predict  near-field  fountain  reingestion  are 
summarized  below. 

1)  Separation  of  the  spreading  ground  jet  is  predicted 
using  Abbott's  data  correlation. 

2)  The  dividing  streamline  in  the  ground  plane  is  predicted 
by  potential  flow. 

3)  The  exhaust  cloud  upper  surface  is  predicted  by  Abbott's 
and  Colin's  data. 

4)  The  exhaust  cloud  is  checked  for  intersection  with  inlet 
streamtubes . 

2.2.3  TR  and  TV  System  Performance  Program 

The  TR  and  TV  System  Performance  Program  consists  of  four 
modules  containing  empirical  data  correlations  and  some 
analyses  developed  during  Tasks  1.1  and  1.3.  Inputs  and 
outputs  of  the  program  are  shown  in  Figure  95.  Detailed 
description  of  the  four  performance  modules  follows. 

Internal  Performance Module 


Data  correlations  developed  during  Tasks  1.1  and  1.3  were 
incorporated  into  the  Internal  Performance  Module,  which 
contains  data  correlations  for  the  following  types  of  cruise 
nozzles,  thrust  reversers,  and  vectoring  nozzles. 
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Figure  92:  HEIGHT  OF  FAR  FIELD  FOUNTAIN 
EXHAUST  CLOUD 
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INPUTS: 


1.  TR  or  TV  TYPE  AND  GEOMETRY 

2.  GEOMETRIC  DISCHARGE  ANGLE 
1  NOZZLE  PRESSURE  RATIO 


TR  AND  TV  SYSTEM  PERFORMANCE  PROGRAM 


OUTPUTS: 

1.  TR  AND  TV  INTERNAL  PERFORMANCE  CONSISTING  OF 

Cy,  Cq,  1\  flg,  7?yg,$>0 

Z  ENGINE  STABILITY  MARGIN 


Figure  95'  TR  AND  TV  SYSTEM  PERFORMANCE  PROGRAM 
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1)  Cruise  nozzles 

a)  Conical 

b)  Annular 

c)  Noncircular 

2)  Thrust  reversers 


a) 

Target  (clamshell  and 

annular) 

b) 

Blocker  deflector  and 

blocker  cascade 

Thrust  vectoring  nozzles 

a) 

Single  bearing 

b) 

Three  bearing 

c) 

Spherical  eyeball 

d) 

Lobstertail 

e) 

External  deflector 

4)  Cascade  lattices  applicable  to  cascade  TR  and  TV 
nozzles 

Data  correlations  for  the  above  nozzles  are  discussed  in 
Sections  2.1  and  2.3.  Computer  program  usage  is  described  in 
Section  3.3. 

Aerodynamic  Interference  Module 

Existing  aerodynamic  data  for  thrust  reverser  systems  are 
for  low-wing,  four-engine  jet  transport  aircraft.  Thrust 
reverser  operation  frequently  caused  unfavorable  aerodynamic 
interference  problems,  such  as  bouyancy  lift  forces  on  the 
wing  and  a  loss  in  the  flap  drag.  Both  effects  were  experienced 
by  the  original  737  clamshell  reverser  deployed  in  front  of 
the  flaps.  When  the  new  target  reverser  was  installed  behind 
the  wing  flaps,  favorable  interference,  i.e.,  decreased  lift, 
increased  drag,  resulted. 

High-wing  transport  aircraft  with  pylon-mounted  engines 
experience  an  increase  in  flap  and  spoiler  drag  due  to  reverser 
operation,  as  shown  in  Figure  96.  The  data  correlation  from 
Ref.  44  gives  the  percent  reduction  in  drag  as  a  function 
of  reverser  momentum  coefficient. 

Thrust  reversers  for  STOL  transports  will  be  designed  to 
deflect  the  flow  forward  and  up  to  avoid  aerodynamic  inter¬ 
ference  and  reingestion  problems.  Existing  data  are  inadequate 
to  predict  aerodynamic  interference  for  configurations  of 
this  type.  Data  are  required  showing  the  effects  of  TR 
type,  location,  forward  speed,  pressure  ratio,  and  discharge 
pattern. 
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PERCENT  AIRPLANE  DRAG  LOSS 
DUE  TO  REVERSE  THRUST  OPERATIONS 


AIRPLANE  REVERSED  AIR  MOMENTUM  COEFFICIENT  •  £  (MV  +  A  PA)R  /q^S* 


Figure  96:  EFFECT  Of  REVERSE  THRUST  ON  AIRPLANE 
DRAG  DURING  GROUND  ROLL 
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Logic  has  been  formulated  in  TEM-357  that  allows  easy 
incorporation  of  aerodynamic  data  into  the  Aerodynamic 
Interference  Module  as  data  become  available. 

*  Reingestion  Module 

The  purpose  of  the  Reingestion  Module  is  to  predict  inlet 
temperature  rise  and  total  pressure  distortion  as  a  function 
of  thrust  reverser  type,  location,  forward  speed,  and  dis¬ 
charge  pattern.  A  vast  amount  of  reingestion  data  have 
been  accumulated  for  commercial  transport  thrust  reversers. 
However,  the  data  are  of  little  value  to  candidate  STOL 
transport  thrust  reverser  installations  because  of  significant 
differences  from  commercial  transport  TR  designs. 

,  Logic  has  been  formulated  in  TEM-357  that  allows  easy  incorpor¬ 

ation  of  reingestion  data  into  the  Reingestion  Module  as  data 
become  available. 

Engine  Stability  Margin  Module 

I  The  purpose  of  the  Engine  Stability  Margin  Module  is  to 

predict  effects  of  thrust  reverser  area  mismatch,  engine 
transient  response,  and  inlet  pressure  and  temperature 
distortion  on  engine  stability  margin.  The  engine  stability 
margin  analysis  was  performed  by  Pratt  &  Whitney  Aircraft. 

The  results  have  been  tabulated  and  incorporated  into  the 

j  Engine  Stability  Margin  Module.  The  analysis  is  summarized 

in  the  following  paragraphs  and  described  in  detail  in 
Ref.  45.  A  sample  case  is  provided  in  Appendix  II. 

Stability  margin  studies  were  performed  on  three  Pratt  & 

Whitney  Aircraft  turbofan  engines  u^ing  a  transient  computer 

I  program.  These  engines  were  extracted  from  the  Pratt  &  Whitney 

Aircraft  file  of  study  engine  cycles  and  consist  of  a  2.0  BPR 
mixed  flow,  and  a  6.0  and  12.0  nonmixed  flow  design,  each 
with  a  fanhigh  compression  system.  These  cycles  are  typical 
of  those  considered  at  different  BPR  levels  for  STOL 
application.  The  objective  of  these  stability  margin  studies 

I  was  to  obtain  a  set  of  data  which  would  permit  evaluation  of 

engine  stability  margin  during  normal  STOL  operation  as  a 
function  of  distortion,  reverser  operation  and  thrust  level 
to  ultimately  define  the  distortion  tolerance  for  each  engine. 
The  reverser  effect  on  the  engine  was  defined  to  be  a  change 
in  exhaust  nozzle  effective  area. 

The  assumed  exhaust  nozzle  area  variation  is  shown  in  Figure 
97.  The  definition  of  the  deployment  characteristic  was 
selected  because  a  review  of  several  thrust  reversers 
showed  an  over  area  of  at  least  20  percent  sometime 
during  deployment.  The  review  also  showed  that  an  under  area 
usually  occurs  near  the  end  of  the  deployment  stroke.  The 
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Figure  97:  AIRFLOW  MA  TCH  CUR  VE  REPRESENTING 
REVERSER  DEPLOYMENT 
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deployment  rate  in  all  cases  was  one  second  from  stowed 
to  fully  deployed.  This  rate  was  identified  as  a  goal  for 
thrust  reverser  design  studies  conducted  in  Part  IB  of 
the  program. 

Data  were  calculated  for  both  steady  and  transient  operation 
at  a  sea  level  flight  condition.  Steady  state  data  show  the 
effect  of  thrust  reverser  area  variation  on  the  stability  margin 
of  the  fan  and  compressor,  as  well  as  the  thrust  split  between 
the  flowstream.  The  thrust  split  can  be  used  to  estimate  the 
reverse  thrust  of  different  types  of  reversers.  The  steady 
state  data  were  obtained  from  each  engine  model  by  simulating 
the  effect  of  a  thrust  reverser  as  a  change  in  either  primary 
or  duct  stream  effective  area.  The  effective  area  values 
investigated  were  the  extremes  of  minimum  and  maximum  area 
positions  as  identified  by  Boeing. 

Transient  operation  consisted  of  engine  accels  and  decels  at 
various  area  values,  as  well  as  thrust  reverser  transitions. 
Accels  and  decels  were  obtained  to  determine  the  influence 
of  off  design  area  value  on  stability  margin  as  well  as  the 
time  it  takes  to  achieve  a  desired  thrust  change.  Thrust 
reverser  transitions  were  obtained  for  the  one  second  deploy¬ 
ment  rate  with  engines  set  at  80  percent  of  takeoff  power, 
and  were  made  from  stowed  to  deployed  as  well  as  from  deployed 
to  stowed. 

It  is  concluded  from  analysis  of  these  data  that  the  influence 
of  a  thrust  reverser  on  component  stability  margin  during  trans¬ 
ient  operation  is  approximately  equivalent  to  changes  obtained 
during  steady  operation  at  the  same  nozzle  positions.  This 
can  be  observed  by  examination  of  accel  and  decel  data,  as 
well  as  thrust  reverser  deployment  data.  This  conclusion 
simplifies  thrust  reverser  analysis  because  it  eliminates  the 
need  for  incorporating  transient  stability  margin  data 
into  the  analysis. 

A  further  simplification  was  derived  for  accel  and  decel 
influence  on  high  compressor  stability  margin.  These  data 
have  been  reduced  to  a  set  of  curves  that  define  the 
loss  or  gain  in  stability  margin  as  a  function  of  engine 
power . 

Other  factors  that  influence  stability  margin  are  pressure  and 
temperature  distortion.  Data  were  obtained  for  each  engine 
on  a  "change  in  stability  margin"  basis.  The  effect  of  dis¬ 
tortion  was  obtained  using  parallel  compressor  theory  re¬ 
sulting  from  180  degree  square  wave  patterns. 


143 


2.3  Task  1.3--Plan  and  Conduct  Supplemental  Tests 

Task  1.3  consisted  of  planning  and  conducting  supplemental 
tests  in  conjunction  with  development  of  computer  programs 
for  thrust  reverser  and  thrust  vectoring  systems.  The  purpose 
of  the  tests  was  to  fill  data  voids  discovered  in  the  liter¬ 
ature  that  would  have  impeded  development  of  computer  programs. 

2.3.1  Identification  of  Technology  Voids 

Boeing  considered  identification  of  data  voids  to  be  an 
important  task  of  this  program.  Computer  programs  developed 
during  Task  1.2  are  based  on  results  of  the  literature  review. 
Also ,  voids  discovered  formed  the  basis  for  supplemental 
tests.  Consequently,  considerable  effort  was  made  in  re¬ 
viewing  data  to  accurately  identify  data  voids  and  justify 
selection  of  models  for  supplemental  tests. 

A  guideline  established  early  in  the  planning  of  supplemental 
tests  was  that  the  tests  would  supply  data  required  to  complete 
development  of  the  Internal  Performance  Module  of  the  TR  and 
TV  System  Performance  Program  (sec  Section  2.2.3).  Obviously, 
data  voids  exist  for  new  thrust  reverser  and  thrust  vectoring 
concepts.  However,  it  is  believed  that  model  testing  of  new 
concepts  should  be  deferred  to  Part  1C  after  adequate  design 
and  evaluations  have  been  completed. 

The  data  review  resulted  in  identification  of  tae  following 
data  voids  in  the  open  literature: 

Blocker  Door  Geometry 

There  was  a  lack  of  a  consistent  set  of  data  that  isolate 
effects  of  blocker  door  geometry  on  the  performance  of  blocker/ 
deflector  and  blocker/cascade  thrust  reversers. 

Multibearing  Vectoring  Nozzles 

There  was  a  lack  of  data  for  multibearing  vectoring  nozzles 
that  show  the  effect  of: 

o  Duct  contraction  ratio 
o  Turning  radius 

o  Low  nozzle  pressure  ratios  (*1.5) 

Also,  there  was  no  data  that  evaluate  the  multibearing 
vectoring  nozzle  as  a  combined  thrust  reverser  and  vectoring 
system  (120  degree  deflection) .  Available  data  for  other  types 
of  thrust  reverser  and  vectoring  systems  are  adequate  to 
formulate  the  Internal  Performance  Module  of  the  TR  and 
TV  System  Performance  Program. 
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2.3.2  Supplemental  Static  Tests 

On  the  basis  of  the  data  review,  a  test  plan  (Ref.  46)  was 
formulated  for  supplemental  static  tests  to  determine: 

1)  Multibearing  vectoring  nozzle  performance  as  a 
function  of  parametric  geometry  variations. 

2)  Blocker  door  geometry  effects  on  performance  of  blocked/ 
cascade  and  blocker A!>f lector  thrust  reversers. 

The  tast  plan  was  submitted  to  the  Air  Force  program  manager 
in  September  1971  and  was  approved  in  October.  Static  test¬ 
ing  was  conducted  from  20  December  1971  to  14  February  1972 
on  the  thrust  vector  rig  at  the  Boeing  Propulsion/Noise 
Laboratories . 

A  test  report  (Ref.  47)  was  prepared  and  submitted  to  the 
Air  Force  program  manager  in  March  1972.  Descriptions  of 
test  models  and  results  are  summarized  in  the  following 
sections.  The  test  results  have  been  incorporated  into  the 
Internal  Performance  Module  of  ?R  and  TV  Systerm  Performance 
Program  TEM-357. 

Multibearing  Vectoring  Nozzle  Static  Test 

The  objectives  of  the  multibearing  vectoring  nozzle  static 
tests  were  to: 

1)  Evaluate  the  vectoring  nozzle  internal  performance  (thrust 
and  mass  flow)  for  variations  of: 

a)  Bearing  angle 

b)  Duct  turning  Mach  number  (inlet  and  exit) 

c)  Duct  turning  radius 

d)  Overall  duct  and  nozzle  length 

e)  Thrust  vector  angle 

f)  Nozzle  pressure  ratio 

2)  Evaluate  the  multibearing  vectoring  nozzle  as  a  thrust 
reverser  system. 

3)  Establish  design  criteria  for  full-scale  nozzle  design. 

A  total  of  10  multibearing  nozzle  test  configurations  were 
tested  including  the  following  geometric  variations: 
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1)  Bearing  angle,  0=  65°,  60° 

2)  Duct  contraction  ratio: 

0  =  65°,  Ajj/Ag  *  1.345,  1.760,  2.94 
0«  60°,  Ad/A£  *  1.345,  2.94 

3)  Duct  turning  radius  ratio: 

0=  65°,  (r/D)min  *  0.641,  0.926 
0=  60°,  (r/D)min  «  0.645,  0.840 

4)  Overall  duct  and  nozzle  length: 

0s  65°,  L/D  =  2.415,  1.996 

0=  60°,  L/D  =  2.944,  2.301 

5)  Thrust  vector  angle 

0  *  65°,  <T=  0°,  30°,  60°  100° 

0  s  0°,  30°,  60* ,  90°,  120° 


In  addition,  one- of  the  models  tested  has  a  converging  duct 
section  that  accelerated  flow  through  the  duct  turn.  The 
model  had  a  contraction  ratio  of  A^/'A^  *  2.94  and  bearing 

plane  angles  of  60  degrees.  All  models  were  tested  over 
the  pressure  ratio  range  from  1.1  to  3.0.  A  test  model  is 
shown  installed  on  the  vectoring  rig  in  Figure  98. 

Many  of the_ configurations  ha$  nearly  equivalent  performance. 
Therefore  to  improve  readability,  data  points  are  not  shown 
on  the  following  comparison  plots.  Plotted  data  for 
each  nozzle  configuration  are  presented  in  Appendix  I  of 
Ref.  47.  Results  of  the  multibearing  static  test  are 
summarized  below. 

1)  Vector  efficiency  and  airflow  match  are  decreased  as 
duct  contraction  ratio  decreases  (and  duct  Mach  number 
increases)  especially  below  contraction  ratio  cf  1.760. 
Airflow  match  was  significantly  lower  for  contraction 
ratio  of  1.345  at  maximum  deflection  angle  (Figure  99 
through  101). 

2)  The  duct  turning  radius  was  not  a  significant  influence 
on  vectoring  nozzle  performance.  Generally,  vectoring 
efficiency  was  improved  0.5  to  1.5  percent  and  airflow 
match  was  improved  by  two  percent  when  the  turning  radius 
ratio  was  increased  at  maximum  deflection  (Figure  102 
and  103) . 
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Figure  98 :  MUL  T! BEARING  VECTORING  NOZZLE  MODEL  INSTALLA TION 
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Figure  99:  (Concluded)  EFFECT  OF  DUCT  CONTRACTION  RA  TIO  ON 

VECTOR  EFFICIENCY  d>  =  65° 
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Figure  101:  EFFl'CT  OF  DUCT  CONTRACTION  RATIO  AT Pj/P^  1.60 
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Figure  133:  1 Concluded, )  EFFECT  OF  DUCT  TURNING  RADIUS  ON 

AIRFLOW  MA  TCH.  6  =  65° 


3)  The  effect  of  overall  duct  length  Indicated  duct  and 
nozzle  length  can  be  made  as  short  as  possible  without 
severe  penalty  to  vectoring  nozzle  performance  (Figures 
104  and  105). 

4}  Reverser  efficiency  of  47  to  49  percent  was  obtained  when 
the  60  degree  bearing  plane  angle  configurations  were 
tested  at  120  degree  deflection  angle.  The  data  indicate 
that  satisfactory  levels  of  reverser  performance  are 
possible  when  the  multibearing  nozzle  is  used  as  a 
thrust  reverser  (Figure  106). 

5)  The  converging  duct  model  configuration  had  equivalent 
vectoring  performance  when  compared  to  a  constant  area 
duct  configuration  with  the  same  contraction  ratio.  This 
design  is  attractive  because  it  is  shorter  than  constant 
area  duct  designs  and  therefore  will  have  less  weight 
(Figuresl07  and  108). 

Thrust  Reverser  Blocker  Door  Geometry  Static  Test 

The  purpose  of  this  static  test  was  to  investigate  blocker 
door  geometry  effects  on  thrust  reverser  performance.  Test 
objectives  were: 

1)  Determine  the  effect  of  blocker  door  deflection  angle 
and  setback  distance  on  performance  of  blocker/deflector 
and  blocker/cascade  thrust  reversers  installed  in  the 
following  ducts: 

a)  Annular  duct  of  a  fan  nozzle 

b)  Circular  duct  of  a  primary  (or  mixed  flow  nozzle) 

2)  Establish  design  criteria  for  the  above  thrust  reverter 
installations. 

Initial  planning  for  this  test  included  four  types  of  thrust 
reversers: 

o  Annular  duct  blocker/cascade 
o  Annular  duct  blocker/deflector 
o  Circular  duct  blocker/cascade 
o  Circular  duct  blocker/deflector 

However,  static  testing  of  the  blocker  door  geometry  model  was 
terminated  after  completion  of  the  annular  duct  blocker/cascade 
model  testing.  It  was  decided  that  the  results  obtained  would 
satisfy  the  ovejrall  objectives  of  the  program  without  expendi¬ 
ture  of  funds  bisyond  the  budgeted  amount.  Since  effects  of 
blocker  door  geometry  were  found  to  be  relatively  small,  lack 
of  data  for  untested  configurations  will  not  significantly 
affect  the  accuracy  of  internal  performance  predictions  of 
the  TR  and  TV  System  Performance  Program. 
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Figure  106:  EVALUATION  OF  THRUST  REVERSER 
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107:  EFFECT  OF  ACCELERA  TING  DUCT  FLOW 
ON  VECTOR  EFFICIENCY 


oO€  *  o 


The  blocker  door  geometry  static  test  model  consisted  of  an 
annular  flow  duct  simulating  a  fern  exhaust  duct,  a  simulated 
engine  strut  (internal),  and  a  cascade  reverser  model  similar 
to  the  type  used  on  the  fan  reverser  on  the  Model  747  JT9D-3A 
engine.  A  photograph  of  the  model  installation  is  shown  in 
Figure  109. 


Inserts  for  45,  90,  and  135  degrees  were  used  to  simulate  the 
blocker  door  geometry.  The  model  was  tested  with  setback 
distances  of  s/h  =  0  and  0.50. 

The  effects  of  blocker  door  geometry  were  evaluated  as  incre¬ 
ments  relative  to  the  0  =  90  degrees,  s/h  =  0,  blocker  door 
configuration.  Data  points  are  not  shown  in  the  comparison 
curves.  Test  data  for  each  model  configuration  showing 
absolute  levels  are  shown  in  Appendix  I  of  Ref.  47.  The  re¬ 
sults  of  the  test  are  summarized  as  follows: 

1)  Data  indicate  optimum  blocker  door  design  for  an  annular 
cascade  reverser  is  the  90  degree  door  geometry.  Reverser 
performance  was  lower  for  the  45  and  135  degree  blocker 
door  geometry  configurations.  (Figure  110) 

2)  Blocker  door  setback  distance  has  minimal  effect  on 
reverser  performance,  as  shown  in  Figure  111.  The  data 
show  that  small  positive  or  negative  changes  in  reverser 
efficiency  occurred  when  the  blocker  door  was  set  back 
from  the  cascade  depending  on  the  nozzle  pressure  ratio. 

Details  of  test  models  and  installation,  instrumentation,  pro¬ 
cedure,  and  results  of  the  multibearing  vectoring  nozzle  and 
blocker  door  geometry  static  tests  are  presented  in  Ref.  47. 
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Figure  111:  EFFECT  OF  BLOCKER  DOOR  SETBACK 
DISTANCE  ON  CASCADE  REVERSER 
PERFORMANCE 
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SECTION  III 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  tasks  conducted  during  Part  1A  of  this  program  have  pro¬ 
duced  several  useful  tools  for  thrust  reverser  and  vectoring 
nozzle  design  studies: 

1)  Thrust  Reverser  and  Thrust  Vectoring  Literature  Review 
Document  (Ref.  1  )  — to  allow  a  reader  to  easily  determine 
sources  of  data  related  to  his  particular  interest. 

2)  Jet  Trajectory  and  Spreading  Program — to  predict  the 
trajectory  and  shape  of  the  thrust  reverser  or  vectoring 
nozzle  exhaust  plume. 

3)  Reingestion  Prediction  Program — to  predict  the  on-set 

of  reingestion  for  arbitrary  thrust  reverser  and  airplane 
configurations  as  a  function  of  geometry  and  flow  con¬ 
dition. 

4)  TR  and  TV  System  Performance  Program — to  predict  TR  and 
TV  nozzle  performance  and  engine  stability  margin 
characteristics . 

In  addition,  static  tests  were  conducted  to  determine: 

1)  Multibearing  vectoring  nozzle  performance  as  a  function 
of  parametric  geometry  variations. 

2)  Blocker  door  geometry  eifects  on  the  performance  of 
annular  blocker/cascade  thrust  reversers. 

The  test  results  filled  data  voids  discovered  during  the 
literature  review  and  were  incoporated  into  the  TR  and  TV 
System  Performance  Program.  The  primary  conclusions  to 
be  drawn  from  the  Part  1A  analytical  work  are: 


1)  The  analytical  and  semi-empirical  methods  developed 

provide  relatively  simple  design  tools  to  evaluate  the 
performance  of  many  types  of  thrust  reverser  and  thrust 
vectoring  systems  and  to  determine  potential  exhaust 
flow  interference  and  reingestion  problems.  They  are 
limited  only  by  the  current-state-of-the-art  of  analysis 
and  existing  empirical  data. 
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2)  There  are  many  sources  of  thrust  reverser  reingestion 
and  aerodynamic  interference  data  for  current  CTOL 
transports.  However,  application  of  the  data  to  STOL 
transports  is  limited  because  the  specific  requirements 
of  STOL  transports  cause  the  thrust  reverser  installation 
to  be  highly  configuration  dependent.  The  requirements 
of  exhaust  flow  directional  control  to  minimize  reinges¬ 
tion  result  in  thrust  reverser  types  that  are  significant¬ 
ly  different  than  conventional  designs.  The  aircraft 
high  lift  system  has  a  direct  influence  on  the  nacelle 
location,  reverser  type,  and  hence  reingestion  and 
aerodynamic  interference  characteristics. 

3)  The  methods  used  to  predict  the  onset  of  reingestion  would 
benefit  from  comparisons  with  experimental  results  that 
isolate  the  specific  type  of  reingestion  i.e.,  crossflow, 
self  (closed  loop),  near  field  fountain,  or  far  field 
fountain. 

4)  The  influence  of  thrust  reverser  operation  on  engine 
stability  margin  during  transients  is  approximately 
equivalent  to  changes  obtained  during  steady  state 
operation  at  the  corresponding  partially  deployed  reverser 
position. 

5)  Existing  data  for  flat  plate  and  curved  external  deflectors 
is  inadequate  to  predict  effects  of  setback  distance  and 
door  length  for  other  types  of  deflector  nozzles.  Static 
tests  were  conducted  during  Part  1C  on  a  hinged  deflector 
nozzle  and  data  correlation  showing  effects  «f  setback 
distance  were  incorporated  into  the  TR  and  TV  System 
Performance  Program. 


Recommendations  for  future  thrust  reverser  and  thrust  vector¬ 
ing  work  pertinent  to  the  tools  developed  for  Part  1A  include 
the  following: 

1)  Interactions  between  thrust  reverser  and  inlet  flow  fields 
and  the  effects  on  airplane  aerodynamic  characteristics 
should  be  studied  using  three-dimensional  potential  flow 
techniques. 

2)  The  closed  loop  reingestion  criteria  should  be  improved 
by  using  a  free  streamline  analysis  or  by  experiment. 

A  parametric  test  of  variables  influencing  self  re¬ 
ingestion  is  recommended. 

3)  Improvements  should  be  made  to  the  Jet  Trajectory  and 
Spreading  Program  to  handle  multiple,  merging  jets,  includ 
ing  effects  of  inlet  suction  and  nonuniform  flow  field. 
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A  low  speed  wind  tunnel  test  should  be  conducted  to 
obtain  thrust  reverser  and  thrust  vectoring  performance 
data  allowing  correlations  of  aerodynamic  interference 
and  reingestion  as  a  function  of  TR  and  TV  system  type, 
location,  forward  speed,  pressure  ratio,  vector  angle, 
and  angle  of  attack  and  yaw. 


SECTION  IV 


COMPUTER  PROGRAM  USAGE 


Three  computer  programs  have  been  developed: 

1)  Jet  Trajectory  and  Spreading  Program  TEM-356A 

2)  Reingestion  Prediction  Program  TEM-356B 

3)  TR  and  TV  System  Performance  Program  TEM-357 

The  programs  are  coded  in  FORTRAN  IV  language  for  the  Control 
Data  Corporation  6600  (131K)  digital  computer.  They  may  be 
compiled  with  either  FORTRAN  IV  (Extended)  or  FORTRAN  IV  (Run) 
compilers  by  changing  the  end-of-file  test  in  input  routines. 
The  standard  FORTRAN  system  library  (FTN)  is  used.  There  is 
no  overlaying  or  use  of  scratch  tapes.  Control  of  the  com¬ 
puter  during  checkout  was  by  the  KRONOS  operating  system. 

Card  deck  arrangements  for  the  programs  are  shown  in  Figures 
112  and  113.  Each  program  consists  of  control  cards,  either 
a  symbolic  or  binary  deck,  and  a  group  of  data  cards.  The 
control  cards  shown  are  common  for  the  three  programs  but  are 
characteristic  of  the  particular  computer  installation.  Data 
cards  for  the  particular  problem  and  the  first  two  control 
cards,  which  contain  run  priority,  time  estimate,  and  user 
identification  are  inserted  into  the  deck.  The  other  cards 
normally  remain  unchanged. 

The  remainder  of  this  section  describes  mechanics  of  input 
data  preparation.  Complete  descriptions  of  card  inputs  to 
each  program  are  presented.  Appendix  II  contains  sample 
cases  to  acquaint  users  with  input  card  format  and  printout 
format  for  each  program.  These  sample  cases  are  also  in¬ 
tended  for  use  in  checking  out  the  programs  on  computer 
facilities. 

3.1  Jet  Trajectory  and  Spreading  Computer  Program  Usage 

Program  TEM-356A  provides  a  simple  method  of  obtaining  jet 
plume  shapes  for  four  types  of  initial  jet  cross  ections. 


TYPE  = 

Initial  Jet  Cross  Section 

Page 

1. 

Circular 

174 

?. 

Rectangular 

176 

3. 

Two-dimensional 

178 

4. 

Annular 

180 
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Figure  f 12:  DECK  A  RRANGEMENT  FOR  PROGRAMS  T EM-358  A,  TEM-356B, 
AND  T EM-357  WHEN  USING  A  SYMBOLIC  DECK 
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I 

Figure  1 13:  DECK  ARRANGEMENT  FOR  PROGRAMS  T EM-356 A,  TEM-356B, 
AND  TEM-357  WHEN  USING  A  BINARY  DECK 
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Data  Input  Format 


All  input  to  program  TEM-356A  is  in  the  "NAMELIST"  format, 
except  one  table  of  values  used  for  TYPE  *=  4.  (page  181). 

The  following  rules  must  be  followed  to  use  NAMELIST  format. 

1)  Card  Column  1  must  be  blank,  always. 

2)  A  data  set  begins  with  a  $NEW  card  and  ends  with  a  $END 
card. 

3)  Commas  are  used  to  separate  numbers  and  variable  names. 

4)  Variable  names  must  be  properly  spelled. 

5)  Card  Columns  1  to  80  are  read  using  NAMELIST.  This 
means  that  CC  73-80  cannot  be  used  for  identification 
purposes . 

6)  Variable  names  and  values  are  read  "free  field."  This 
means  imbedded  blanks  may  be  used  at  the  user's  dis¬ 
cretion,  and  that  more  than  one  variable  may  be  input 
per  card.  However ,  the  input  system  described  on  the 
following  pages  is  recommended  because  it  separates 
the  numbers  and  makes  data  checking  easier. 

The  general  card  stacking  arrangement  is  shown  in  Figure  114. 
The  user  may  run  as  many  jet  plumes  as  he  desires  in  a 
single  computer  submission.  There  are  no  restrictions  as  to 
sequence  of  the  four  types.  Detailed  card  inputs  for  the  four 
types  are  contained  in  the  next  section.  Users  are  cautioned 
to  check  their  inputs  carefully  before  submitting  a  computer 
run. 
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Jet  Trajectory  and  Spreading  Program  Card  Input 


TYPE  «  1.  Figure  115  displays  the  data  card  arrangement 
for  circular  jet  plumes.  A  description  of  the  card  input 
follows: 


Card 

Column 

Code 

Explanation 

■ 

$NEW 

m 

$NEW 

$NEW  must  be  punched  in 

CC  2-10. 

Card  1 

2-10 

TYPE  « 

TYPE  *  must  be  punched  in 

CC  2-10. 

11-20 

1., 

indicates  circular  jet 

Card  2 

2-10 

XYZ0  ** 

XYZ0  -  must  be  punched  in 

CC  2-10. 

11-20 

Xo' 

jet  origin  in  reference  co¬ 

u 

ordinate 

21-30 

V 

system 

31-40 

V 

Note:  The  fourth  character 

Q 

in  XYZ0  is  a  "zero",  not  the 
letter  "O". 

Card  3 

2-10 

TXYZ  « 

TXYZ  »  must  be  punched  in 

1 

CC  2-10. 

11-20 

vector  components  defining 

initial 

21-30 

V 

jet  direction 

31-40 

V 

Card  4 

2-10 

D  « 

D  ■  must  be  punched  in  CC  2-10. 

11-20 

d, 

jet  Initial  diameter 

Card  5 

2-10 

OELSD  - 

As/d, 

OELSD  ■  must  be  punched  in 

CC  2-10. 

11-20 

arc  length  spacing/initial 
diameter 

Card  6 

2-10 

Q RAT 10  « 

QRATI0  ■  must  be  punched  in 

CC  2-10. 

11-20 

/«u> 

jet  to  freestream  dynamic 
pressure  ratio 

Card  7 

2-10 

UINF  « 

UINF  ■  must  be  punched  in  CC  2-10. 

11-20 

uoo, 

freestream  velocity 

174 


Figun  1 1S  DA  TA  CARD  ARRANGEMENT  FOR  TYPE  •  1.  CIRCULAR  JET 
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TYPE  =  2.  Figure  116  displays  the  data  card  arrangement  for 
rectangular  jet  plumes.  A  description  of  the  card  input  follows: 


Card 

Column 

Code 

Explanation 

$NEW 

2-10 

$NEW 

$N£W  must  be  punched  in  CC  2-10. 

Card 

1 

2-10 

TYPE  « 

TYPE  *  must  be  punched  in  CC  2-10. 

11-20 

2., 

indicates  rectangular  jet 

Card 

_2 

2-10 

XYZ0  * 

XYZ0  must  be  punched  in  CC  2-10. 

11-20 

xo' 

jet  origin  in  reference  co¬ 

ordinate 

■'IBM 

yo' 

system 

1  I 

zo' 

Note:  The  fourth  character  m 

o 

XY70  is  a  "zero",  not  the 
letter  "O". 

Card 

2-10 

TXYZ  « 

TXYZ  »  must  be  punched  in 

CC  2-10. 

11-20 

V 

vector  components  defining 

X 

initial 

21-30 

V 

jet  direction 

31-40 

V 

Card 

_4 

2-10 

A  » 

A  *  must  be  punched  in  CC  2-10. 

i 

11-20 

i 

a, 

longitudinal  dimension  of 
initial  cross  section 

Card 

_5 

2-10 

B  - 

B  *  must  be  punched  ir.  CC  2-10. 

11-20 

1 

b. 

lateral  dimension  of  initial 
cross  section 

Card 

_6 

2-10 

DELSD  - 

DELSD  ■  must  be  punched  in 

CC  2-10. 

11-20 

As/de, 

arc  length  spacing/equivalent 
diameter 

Card 

2 

2-10 

QRATI0  » 

QRATI0  ■  must  be  punched  in 

CC  2-10. 

11-20 

jet  to  freestream  dynamic  pressure 
ratio 

Card 

2 

2-10 

UINF  ■ 

UINF  ■  must  be  punched  in  CC  2-10. 

11-20 

freestream  velocity 

Card 

2 

2-10 

ALPHA  - 

ALPHA  ■  must  be  punched  in  CC  2-10. 

11-20 

angle  of  attack 
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TYPE  =  3.  Figure  117  displays  the  data  card  arrangement  for 
two-dimensional  jet  plumes.  A  description  of  the  card  input 
follows: 


C>rd 


$NEW 
Card  1 


Card  2 


Card  3 


Card  4 


Card  5 


Card  6 


Card  7 


Card  8 


Card  9 


Column 

Code 

Explanation 

2-10 

$JVEw 

$NEW  must  be  punched  in  CC  2-10. 

2-10 

TYPE  « 

i 

TYPE  ■  must  be  punched  in 

CC  2-10. 

11-20 

3., 

indicates  two-dimensional  cross 
section 

2-10 

XYZ0  - 

XYZ0  -  must  be  punched  in 

CC  2-10. 

11-20 

V 

jet  origin  in  reference  co¬ 
ordinate 

!  21-30 

31-40 

V 

system 

V 

Note:  The  fourth  character  in 
XYZ0  is  a  "zero",  not  the 
letter  "O". 

! 

2-10 

TXYZ  » 

TXYZ  *  must  be  punched  in 

CC  2-10. 

11-20 

V 

vector  components  defining 
initial 

21-30 

V 

jet  direction 

31-40 

V 

2-10 

T  » 

T  »  must  be  punched  in  CC  2-10. 

11-20 

t  f 

jet  initial  thickness 

2-10 

W  - 

W  ■  must  be  punched  in  CC  2-10. 

11-20 

W, 

jet  initial  width 

2-10 

DELSD  - 

DELSD  •  must  be  punched  in 

CC  2-10. 

11-20 

As/t, 

arc  length  spacing/jet  initial 
thickness 

2-10 

Q RAT 10  « 

QRAT10  ■  must  be  punched  in 

CC  2-10. 

11  20 

jet  to  freestream  dynamic 
pressure  ratio 

2-10 

ALPHA  - 

ALPHA  ■  must  be  punched  in 

CC  2-10. 

11-20 

angle  of  attack 

2-10 

PSI  » 

PSI  *  must  be  punched  in  CC  2-10. 

11-20 

>> 

angle  of  yaw 

178 
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Figure  1 17:  DA  TA  CARD  ARRANGEMENT  FOR  TYPE  -  X  TWO  DIMENSIONAL  JET 
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TYPE  =  4.  Figure  118  displays  the  data  card  arrangement 
for  annular  jet  plumes.  A  description  of  the  card  input 
follows: 


Card 


$NEW 


Card  1 


Card  2 


Card  3 


Card  4 


Card  5 


Card  6 


Card  7 


Card  8 


Card  9 


Column 

Code 

Explanation 

2-10 

$NEW 

$N£W  must  be  punched  in  CC  2-10. 

2-10 

TYPE  = 

TYPE  =  must  be  punched  in 

CC  2-10. 

11-20 

4., 

indicat  ;s  annular  jet  cross 
section 

2-10 

XYZA  = 

XYZA  -  must  be  punched  in  CC  2-10. 

Xa' 

origin  of  annular  jet  coordinate 

Cfc 

V 

system 

V 

2-10 

AXYZ  * 

AXYZ  *  must  be  punched  in  CC  2-10. 

11-20 

V 

vector  components  defining 
cruise 

21-30 

V 

no22le  direction 

31-40 

V 

2-10 

T  - 

T  *  must  be  punched  in  CC  2-10. 

11-20 

t, 

jet  initial  thickness 

2-10 

R  - 

R  ■  must  be  punnhed  in  CC  2-10. 

11-20 

r, 

annular  jet  initial  radius 

2-10 

DELSD  » 

DELSD  ■  must  be  punched  in 

CC  2-10. 

11-20 

As/t 

arc  length  spacing/jet  initial 
thickness 

2-10 

MR  - 

MR  »  must  be  punched  in  CC  2-10. 

11-20 

MR, 

number  of  annular  jet  arms 

2-10 

QRATI0  - 

QRATX0  *>  must  be  punched  in 

CC  2-10. 

11-20 

jet  to  frees tream  dynamic 
pressure  ratio 

2-10 

ALPHA  - 

ALPHA  *  must  be  punched  in 

CC  2-10, 

11-20 

angle  of  attack 

Card 

Column 

Code 

Explanation 

Card  10 

2-10 

PSI  = 

PSI  =  must  be  punched  in  CC  2-10. 

11-20 

% 

angle  of  yaw 

Card  11 

2-10 

PC0DE  = 

PC0DE  *  must  be  punched  in 

CC  2-10. 

11-20 

PC0DE , 

=  0 . ,  no  punched  cards  produced 
«  1 . ,  punched  cards  produced 

SEND 

2-10 

SEND 

SEND  must  be  in  CC  2-10. 

The  following  additional  data  must  be  input  for  TYPE  =  4., 

Note  that  NAMELIST  format  is  not  used  for  these  inputs.  Input 
format  is  2F10.0  (two  fields,  each  ten  digits  wide). 


Card 

Column 

Code 

Explanation 

Card 

1-10 

6 

arm  clock  angle 

Set  A 

11-20 

a 

flow  turning  angle 

Note:  There  must  be  MR  of  the  cards,  two  numbers  per  card, 
where  MR  is  the  number  of  annular  jet  arms. 


181 


Figure  1 18:  DA  TA  CARD  ARRANGEMENT  FOR  TYPE  =  4.  ANNULAR  JET 
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3.2  Reingestion  Prediction  Computer  Program  Usage 
Data  Input  Format 

Inputs  describing  jet  plume  geometry  use  the  NAMELIST  format 
and  are  nearly  identical  to  that  described  for  the  Jet  Trajectory 
and  Spreading  Program.  The  only  exceptions  are  noted  below: 

1)  ALPHA  is  not  input.  The  theory  assumes  the  frees tream 
flow  is  parallel  to  the  ground  plane  and  automatically 
sets  Of  b  0.  Yaw  angle  y  is  input  as  usual. 

2)  PC0DE  is  not  input.  The  program  assumes  PC0DE  =  0. 

3)  A  new  parameter,  XC0WL,  is  input  for  all  jet  types.  XC0WL 
is  the  length  of  the  cowl  from  the  inlet  leading  edge  to 
the  front  of  the  thrust  reverser  exit.  XC0WL  is  used 

in  place  of  plate  length  in  the  criterion  for  reattachment 
of  exhaust  flow  to  the  cowl  surface  (see  Figure  76). 

4)  A  new  parameter,  N0I,  is  inpui  for  all  jet  types.  N0I 
is  the  inlet  number  associated  with  a  particular  jet. 

5)  UINF  must  be  input  for  all  type  jets.  In  TEM-356A,  UINF 
is  required  for  TYPE  *  1.  and  2.  only. 

6)  The  inlet  data  are  input  in  10F7.0  format. 

The  general  card  stacking  arrangement  is  shown  in  Figure  119. 

The  user  may  submit  up  to  four  inlets  and  four  jets  (other 
than  TYPE  «  4.)  in  a  single  case.  Only  one  TYPE  ■  4.  j»fe-  may 
be  run  per  case.  However,  as  many  cases  as  are  desired  can 
be  stacked  for  a  single  computer  submission. 
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REINGESTION  PREDICTION  PROGRAM  CARD  INPUT 


Figure  119  displays  the  data  card  arrangement  for  the 
Reingestion  Prediction  Computer  Program.  All  numbers  must 
be  punched  with  a  decimal  point.  The  inlet  card  input  data 
are  input  in  7F10.0  format  rather  than  NAMELIST. 


Card 

Cr"i  umn 

Code 

Explanation 

Card  1 

1-10 

NI 

=  number  of  inlets 

1.  £  NI  *  4. 

11-20 

ZGRND 

height  of  (x#y,zi  reference 
coordinate  system  from  ground 
plane 

Card 

Set  2 

t— * 
i 

o 

Xi 

position  of  inlet  center  in 

*i 

reference  coordinate  system 

Zi 

31-40 

di 

inlet  diameter 

41-50 

vhi 

inlet  hilite  velocity 

Note:  There  must  be  NI  cards 

in  Card  Set  2,  five  numbers 
per  card. 

Card  3 

1-10 

NJ  / 

*  number  of  jets 

1.  &  NJ  4. 

Figure  120  displays  additional  r«AMELIST  data  cards  for  use 
in  program  TEM-356B.  The  cards  are  inserted  in  any  order 
between  the  $NEW  and  $END  cards  with  the  other  NAMELIST 
cards  describing  the  jet  geometry.  A  description  of  the 
additional  input  follows: 


Card 

Column 

Code 

Description 

Card  NI 

2-10 

XC0WL  - 

XC0WL  -  must  be  punched  in 

CC  2-10. 

11-20 

xcowl , 

length  of  cowl  surface 

Card  N2 

2-10 

N0I  - 

N0I  ■  must  be  punched  in 

CC  2-10. 

11-20 

N0I, 

inlet  number  associated  with 
this  jet 

Card  N3 

2-10 

UINF  « 

UINF  *  must  be  punched  in 

CC  2-10. 

11-20 

_ 

freestream  velocity,  ft/sec 

Figure  1 19:  DA  TA  CARD  ARRANGEMENT  FOR  REINGESTION 
PREDICTION  PROGRAM  TEM-356B 
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Figure  120:  ADDITIONAL  NAME  LIST  DA  TA  CARDS  FOR 

REINGESTION  PREDICTION  PROGRAM  TEM-356B 
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3.3  TR  and  TV  System  Performance  Computer  Program  Usage 

Program  TEM-357  consists  of  four  modules  as  noted  below.  Four- 
letter  capitalized  control  words  are  used  to  direct  program 
execution  to  the  appropriate  subroutine. 

Internal  Performance  Module  (INTE) 

Engine  Stability  Margin  Module  (EWGI) 

Aerodynamic  Interference  Module  (AERO) 

Reinge&tion  Module  (REIN) 

Data  Input  Format 

Input  to  TEM-357  falls  into  three  categories:  control  cards, 
title  cards,  and  numeric  input.  As  their  name  implies,  the 
control  cards  control  execution  of  the  program.  The  control 
cards  must  be  punched  in  card  columns  1  to  4.  All  input  data, 
except  title  cards  and  control  cards,  are  punched  in  number 
fields  seven  columns  wide,  with  ten  fields  per  card.  Decimal 
points  should  always  be  punched  for  every  input  number. 

The  general  card  stacking  arrangement  is  shown  in  Figure  121. 
There  are  no  restrictions  as  to  sequence  of  the  four  major 
modules.  The  control  word  END  punched  in  card  columns  1  to  3 
terminates  program  execution.  Detailed  card  inputs  for  the 
four  modules  are  contained  in  the  following  sections. 

Internal  Performance  Module  (INTE) 

The  Internal  Performance  Module  contains  data  correlations  for 
the  following  types  of  cruise  nozzles,  thrust  reversers,  and 
vectoring  nozzles. 


cruise 

nozzles 


thrust 

reversers 


thrust 

vectoring* 

nozzles 


cascade  losses f 


Nozzle 

Control  Word 

Page 

conical 

CONI 

190 

annular 

ANNU 

192 

noncircular 

EQUI 

194 

target  (clamshell 

and  annular) 

TARG 

196 

blocker  deflector 

BLOC 

200 

single  bearing 

SING 

202 

three-bearing 

THRE 

204 

spherical  eyeball 

SPHE 

206 

lobstertail 
external  deflector 

LOBS 

208 

cascade  lattices 

CASC 

210 

187 


(Blank) 


188 


i 


Figure  121:  DA  TA  CARD  ARRANGEMENT  FOR  SEVERAL  CASES,  PROGRAM  TEM-357 
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Detailed  card  inputs  for  each  type  of  nozzle  are  described 
on  the  pages  noted  in  the  foregoing  table.  Several  of  the 
subroutines  (CONI,  ANNU,  EQUX,  SING,  THRE,  SPHE,  LOBS)  were 
modified  during  Task  3.2  in  order  to  make  Reynolds  number 
scale  corrections.  The  full  scale  Reynolus  number  (calculated 
at  a  reference  nozzle  pressure  ratio  of  1.5)  is  punched  in  the 
eighth  field  (card  column  50-56)  in  the  E  format  and  must  be 
right  adjusted.  For  example,  a  full  scale  Re  of  10.4x10* 
would  be  punched  as  follows: 


Card  Column 

A. 

S- . 


50 

51 

52 

53 

54 

55 

56 

1 

0 

• 

4 

E 

6 

The  check  mark  (✓')  indicates  a  blank.  Sample  Case  Number 
Two  in  Appendix  II  (page  259)  illustrates  the  scale  correction 
inputs  and  outputs. 

C0NI  Input 


Figure  122  displays  the  data  card  arrangement  for  conical 
cruise  nozzles.  A  description  of  the  card  input  follows. 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

C0NI 

Control  card  -  contains  the 
word  C0NI 

Card  2 

1-7 

LEVEL 

»  1.  indicates  "Level  1"  per¬ 
formance  predictions 

■  2.  indicates  "Level  2"  per¬ 
formance  predictions 

8-14 

15-21 

Ar/l 

OC 

nozzle  offset/duct  length 
internal  wall  angle  at  nozzle 
trailing  edge  in  degrees 

Omit  0<  if  LEVEL  -  1. 

22-28 

VD2 

nozzle  entrance/exit  diameter 

Omit  Dj/Dj  if  LEVEL  -I. 

50-56 

Ref. 

full  scale  Reynolds  number  at  a 
reference  nozzle  pressure  ratio 
of  1.5.  Punched  in  E  Format 
and  right  adjusted.  See  above  • 
example. 

Omit  if  scale  corrections  are  not 
desired. 

Figure  122:  DA  TA  CAHD  ARRANGEMENT  FOR  CONICAL  CRUISE  NOZZLES 


ANNU  Input 

Figure  123  displays  the  data  card  arrangement  for  annular 
cruise  nozzles.  A  description  of  the  card  input  follows. 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

ANNU 

control  card  -  contains  the  word 
ANNU 

Card  2 

1-7 

LEVEL 

*  1.  indicates  "Level  1"  per- 

formance  predictions 

■  2.  indicates  "Level  2"  per- 

formance  predictions 

8-14 

Ar/l 

nozzle  offset/duct  length 

15-21 

<X 

internal  wall  angle  at  nozzle 
trailing  edge  in  degrees.  Omit 

OC  if  LEVEL  -1. 

22-28 

AD/AE 

duct  entrance/nozzle  exit  area 

U  b 

Omit  Ajj/Ag  if  LEVEL  *  i. 

2P-35 

Ro 

nozzle  exit  outer  radius 

36-42 

Ri 

nozzle  exit  inner  radius 

50-56 

Refs 

full  scale  Reynolds  number  at  a 
reference  nozzle  pressure  ratio 
of  1.5.  Punched  in  E  format  and 

right  adjusted.  See  example  on 
page  190.  Omit  if  scale  corrections 
are  not  desired. 

Card  3 

1-7 

N 

■  number  of  pressure  ratios 
in  Card  Set  4.  1.  *  Ns  20. 

Card 

Set  4 

1-7 

8-14 

<  Vp«. » , 
<VP«  >2 

nozzle  pressure  ratios  at  which 
performance  predictions  are  de¬ 
sired.  There  must  be  N  of  these 

15-21 

<pt/p«*  N 

values,  ten  numbers  per  card. 

22-28 

<vp-»; 

29-35 

(VP~>5 

36-42 

<VP»  )c 

43-49 

50-56 

<VP-»8 

57-63 

(PT/P«  > 

64-70 

(PT/P«  ) 

T  _ « 

_ 
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EQUI  Input 

Figure  124  displays  the  data  card  arrangement  for  noncircular 
shaped  cruise  nozzles.  A  description  of  the  card  input 
follows: 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

EQUI 

control  card  -  contains  the 
word  EQUI 

Card  2 

1-7 

LEVEL 

**  1.  indicates  "Level  1"  per¬ 
formance  predictions 

=  2.  indicates  "Level  2"  per¬ 
formance  predictions 

8-14 

Ar/l 

nozzle  offset/duct  length 

15-21 

£X 

internal  wall  angle  at  nozzle 
trailing  edge  in  degrees 

Omit  o<  if  LEVEL  *  1. 

22-28 

Dl/D2 

nozzle  entrance/exit  diameter 

Omit  0^2  if  LEVEL  «  1. 

53-60 

Refs 

full  scale  Reynolds  number  at 
a  reference  nozzle  pressure 
ratio  of  1.5.  Punched  in  E 

) 

format  and  right  adjusted.  See 
example  on  page  190.  Omit  if 

scale  corrections  are  not  desired. 

Card 

Set  4 

1-7 

8-14 

<Vp«’i 
( VP“  >  2 

nozzle  pressure  ratios  at  which 
performance  predictions  are  de- 

sired.  There  must  be  N  of  these 

15-21 

<  VP« ) 

values ,  ten  numbers  per  card. 

22-28 

(PT/Poo  ) 
(Pt/Poo)5 

29-35 

36-42 

<VP«>6 

43-49 

<VP«>7 

50-56 

(PT/P«, ) ' 

<vp->. 

57-63 

64-70 

<Vr« ) 

1  10 

194 


Figure  124:  DA  TA  CARD  ARRANGEMENT  FOR  IRREGULAR  SHAPED  CRUISE  NOZZLES 
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Figures  125a  and  125b  display  the  data  card  arrangement  for 
annular  and  clamshell  target  thrust  reversers.  A  description 
of  the  card  input  follows: 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

TARG 

control  card  -  contains  the 
word  TARG 

Card  2 

1-7 

TAR 

=  1.  indicates  clamshell  target 
thrust  reverser 

=  2.  indicates  annular  target 
reverser 

Card  3A* 

1-7 

blockage  angle  in  degrees. 

see  Figure  23 

8-14 

e 

door  angle  in  degrees,  see 

Figure  23 

15-21 

C/H 

throat  gap/ annulus  height,  see 
Figure  23 

22-28 

VAc 

reverser  exit  area/cruise 
nozzle  exit  area 

Card  3C** 

1-7 

L/D 

door  length/nozzle  exit  diameter 

’  8-14 

X/D 

setback  distance/nozzle  exit 
diameter 

15-21 

LH/D 

average  lip  height/nozzle  exit 
diameter 

22-28 

A 

sweep  angle  in  degrees,  see 

Figure  13 

29-35 

e 

arc  angle  in  degrees,  see  Figure  13 

36-42 

« 

cone  angle  in  degrees,  see 

Figure  13 

43-49 

0 

bevel  angle  in  degrees,  see 

Figure  13 

Note:  *  Omit  card  3A  for  clamshell  target  thrust  reversers. 

**  Omit  Card  3C  for  annular  target  thrust  reversers. 


Card 


Card  4 


Card 
Set  5 


Co! 

umn 

Code 

1- 

■7 

N 

1- 

7 

(PT/P“o 

>1 

8- 

14 

<VP<* 

>2 

15- 

■21 

(?T/Poo 

) 

*5 

22- 

•28 

(PT/P«, 

>4 

29- 

35 

(PT/Poo 

>5 

36- 

42 

(PT/Poo 

>6 

43- 

49 

(PT/Poo 

>7 

50- 

56 

(PT/Poo 

) 

8 

57- 

63 

{PT/P  00 

) 

9 

64- 

70 

(P, p/Poo 

) 

1 

10 

Explanation 


=  number  of  pressure  ratios  in 
Card  Set  4.  1.  -=  N  20. 


nozzle  pressure  ratios  at  which 
performance  predictions  are  de¬ 
sired.  There  must  be  N  of  these 
values,  ten  numbers  per  card. 
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b)  CLAMSHELL  TARGET  THRUST  REVERSER  INPUTS 


Figure  125:  DA  TA  CARD  A  RRANGEMENT  FOR  TARGET  THRUST  REVERSERS 
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Preceding  page  blank 


Figure  126  displays  the  data  card  arrangement  for  blocker 
deflector  and  blocker  cascade  thrust  reversers.  A  description 
of  the  card  input  follows: 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

BL0C 

control  card  -  contains  the  word 
BL0C 

Card  2 

1-7 

BFLAG 

a  0.  indicates  blocker  deflector 

thrust  reverser 

■  1.  indicates  blocker  cascade 

thrust  reverser 

Card  3 

1-7 

eso 

blocker  door  angle  in  degrees, 
see  page  162 

8-14 

s/h 

blocker  door  setback/annulus 
height,  see  page  162.  Omit  if 
BFLAG  ■  0. 

Card  4 

1-7 

N 

■  number  of  pressure  ratios  in 

Card  Set  4.  1.  £  Ns  20. 

Card 

Set  5 

1-7 

8-14 

'2 

nozzle  pressure  ratios  at  which 
performance  predictions  are  de¬ 
sired.  Ther  must  be  N  of  these 

15-21 

<PT/Poo 

>3 

values,  ten  numbers  per  card. 

22-28 

(PT/P« 

>4 

29-35 

(PT/Po. 

>5 

36-42 

(PT/P~ 

>6 

43-49 

<vp- 

50-56 

(PT/P- 

>8 
^  9 

57-63 

<PT/P~ 

64-70 

<PT/P- 

) 

10 
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Figure  126:  DA  TA  CARD  A RRANGEMENT  FOR  BLOCKER  DEFLECTOR  AND 
BLOCKER  CASCADE  THRUST  RE  VERSE RS 
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SING  Input 

Figure  127  displays  the  data  card  arrangement  for  single 
bearing  vectoring  nozzles,  A  description  of  the  card  input 
follows: 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

SING 

control  card  -  contains  the  word 
SING 

Card  2 

1-7 

<x 

bearing  plane  angle  in  degrees, 
see  Figure  28 

8-14 

(9 

bearing  duct  angle  in  degrees, 
see  Figure  28 

15-21 

<f> 

bearing  rotation  angle  in 
degrees,  see  Figure  28 

22-28 

SI 

flow  turning  angle  measured 
from  cruise  nozzle  centerline 

ir.  degrees.  See  page  42. 

Program  automatically  calculates 

A.  if  it  is  not  input. 

29-35 

Ar/l 

nozzle  cffset/duct  length 

50-56 

Refs 

full  scale  Reynolds  number  at 
a  reference  nozzle  pressure 
ratio  at  1.5.  Punched  in  E 

format  and  right  adjusted  as 
shown  on  page  190.  Omit  if 
scale  corrections  are  not  desired. 

Card  3 

1-7 

N 

=  number  of  pressure  ratios 
m  Card  Set  4.  1. S  N  &  20. 

Card 

Set  4 

1-7 

(PT/Poc  ) 

nozzle  pressure  ratios  at  which 

8-14 

(PT/P~  ) 

T  2 

performance  predictions  are  de¬ 
sired.  There  must  be  N  of  these 

15-21 

(PT/P*>  )  3 

values,  ten  numbers  per  card. 

22-28 

29-35 

(PT/P~  ]  4 
(PT/Poo  ) 

36-42 

{P.r/P-o  )  £ 

43-49 

<pT/p~ » 7 

50-56 

<  VP~  1  6 

57-63 

(PT/Poo  )9 

64-70 

(PT/Po*  ) 

1  10 

I 

I 

l 


202 


< 
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THRE  Input 


Figure  128  displays  the  data  card  arrangement  for  three  bearing 
vectoring  nozzles.  A  description  of  the  card  input  follows: 


Card 


Column  Code 


Explanation 


Card  1 


THRE 


control  card  -  contains  the  word 
THRE 


Card  2 


bearing  plane  angle  in  degrees, 
see  page  145 


8- a  4  Ac/At; 

15-21  C T 


duct  entrance/nozzle  exit  area 

vector  angle  in  degrees,  see 
page  145 


22-28 


length  of  center  section 


29-35 


total  nozzle  length 


36-42 


duct  entrance  diameter 


50-56 


full  scale  Reynolds  number  at 
a  reference  nozzle  pressure 
ratio  of  1.5.  Punched  in  E 
format  and  right  adjusted  as 
shown  on  page  190.  Omit  if  scale 
corrections  are  not  desired. 


Card 
Set  4 


(PT/Poo  )  x 

<pT/poo  >2 

(Pt/Poo>3 

(PT/Poo  )  4 
(pT/poo  )  5 

(PT/poo  )  e 
(pT/poo  >7 

50-56  (pT/poo  ) 

57-63  (PT/Pooi® 

64-70  (pT/poo  ) 


8-14 

15-21 

22-28 

29-35 

36-42 

43-49 


nozzle  pressure  ratios  at  which 
performance  predictions  are  de¬ 
sired.  There  must  be  N  of  these 
values,  ten  numbers  per  card. 


SPHE  Input 


Figure  129  displays  the  data  card  arrangement  for  spherical 
eyeball  vectoring  nozzles.  A  description  of  the  card  input 
follows : 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

SPHE 

control  card  -  contains  the  word 

SPHE 

Card  2 

1-7 

0 

mechanical  vector  angle  in  degrees 

50-56 

Refs 

full  scale  Reynolds  number  at  a 

reference  nozzle  pressure  ratio 

of  1.5.  Punched  in  E  format 

and  right  adjusted  as  shown  on 

page  190.  Omit  if  scale  corrections 

are  not  desired. 

Card  3 

1-7 

N 

=>  number  of  pressure  ratios  in 

Card  Set  4.  1.  afc  N  s£.  20. 

Card 

Set  4 

1-7 

(PT/Poo  ) 

nozzle  pressure  ratios  at  which 

Q  «.  1  A 

1  1 
(V  /P  \ 

performance  predictions  are  de- 

O 

*■  CX>  9 

sired.  There  must  be  N  of  these 

15-21 

(PT/Poo  ) 

values,  ten  numbers  per  card. 

22-28 

<VP->4 

29-35 

(PT/Poo  ) 

36-42 

(PT/Poo)^ 

43-49 

(Pt/Ppo)? 

50-56 

<  VP~  )  8 

57-63 

(PT/Po.  )  9 

64-70 

(PT/Poo  ) 

10 

Figure  129:  DATA  CARD  ARRANGEMENT  FOR  SPHERICAL  EYEBALL  NOZZLES 
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Figure  130  displays  the  data  card  arrangement  for  lobstertail 
or  aft-hood  deflector  nozzles.  A  description  of  the  card  input 
follows : 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

L0BS 

control  card  -  contains  the  word 
L0BS 

Card  2 

1-7 

e 

mechanical  vector  angle  in 
degrees 

50-56 

Refs 

full  scale  Reynolds  number  at 
a  reference  nozzle  pressure 
ratio  of  1.5.  Punched  in  E 

format  and  right  adjusted  as 
shown  on  page  190.  Omit  if 
scale  corrections  are  not  de¬ 
sired  , 

Card  3 

1-7 

N 

=  number  of  pressure  ratios 
in  Card  Set  4.  1.  £  N  ^  20. 

Card 

Set  4 

1-7 

8-14 

(PT/Poo  ) 

(PT/POo  )  2 

nozzle  pressure  ratios  at  which 
performance  predictions  are  de¬ 
sired..  There  must  be  N  of  these 

15-21 

(VP“>3 

<vp“>4 

values,  ten  numbers  per  card. 

22-28 

29-35 

<  V1’"’  1  5 

<VP~>6 

3C-42 

43-49 

<vp~>7 

50-56 

(pT/po.)8 

<VP~)9 

57-63 

64-70 

(PT/Poo  ) 

T  10 
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Figure  130: 


DATA  CARD  ARRANGEMENT  FOR  LOBSTERTAIL 
OR  AFT-HOOD  DEFLECTOR  NOZZLES 
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EXTE  Input 


Figure  131  displays  the  data  ,'ard  arrangement  for  external 
deflector  nozzles.  A  description  of  the  card  input  follows: 


Column 


15-21 

22-28 

29-35 

36-42 

43-49 

50-56 

57-63 


Code 


EXTE 


CFLAG 


D  /D 
n 


<1 VP<*  >  x 

<pT/p~  >  2 
(PT/poo  )  3 

(PT/P«  ]  4 

(PT/Poe  }  5 
<PT/Poo)g 

(PT/Poo)^ 

(PT/p<»)8 
{PT/Peo  ) 


Explanation 


control  card  -  contains  the  word 
EXTE 

=  0.  indicates  flat  plate 
deflector 

»  1.  indicates  curved  deflector 

=  2.  indicates  hinged  external 
deflector 

mechanical  vector  angle  in  degrees, 
Omit  if  CFLAG  ■  2. 

hinged  deflector  setback  dis¬ 
tance/nozzle  diameter  in  TV 
mode,  see  page  266.  Omit  if 
CFLAG  =  0.  or  1. 

nozzle  diameter  in  TV  mode/cruise 
nozzle  diameter,  see  page  266. 

Omit  if  CFLAG  =  0.  or  1. 

deflection  angle  defined  in 
sketch  on  page  266.  Omit  if 
CFLAG  *  0.  or  1. 

=  number  of  pressure  ratios  in 
Card  Set  4.  1.  ^  N  &  20. 


nozzle  pressure  ratios  at  which 
performance  predictions  are  de¬ 
sired.  There  must  be  N  of  these 
values,  ten  numbers  per  card. 


63-70  (PT/Poc) 


Figure  131:  DA  TA  CARD  ARRANGEMENT  FOR  EXTERNAL  DEFLECTOR  NOZZLES 
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Figure  132  displays  the  data  card  arrangement  for  cascade  loss 
predictions.  A  description  of  the  card  input  follows: 


Card 

Column 

Code 

-----  . . 

Explanation 

Card  1 

1-4 

CASC 

control  card  -  contains  the  word 
CASC 

Card  2 

1-7 

a. 

inlet  blade  angle  in  degrees, 
see  Figure  38 

8-14 

gas  inlet  angle  in  degrees, 
see  Figure  38 

15-21 

c/s 

blade  chord/pitch  ratio,  see 

Figure  38 

22-28 

o/s 

opening/pitch  ratio,  see  Figure  38 

29-35 

s/e 

pitch/curvature  ratio,  see 

Figure  38 

36-42 

t/c 

thickness/chord  ratio,  see 

Figure  38 

43-49 

Vs 

trailing  edge  thickness/pitch 
ratio,  see  Figure  38 

50-56  . 

T 

tn 

nozzle  total  temperature,  °R 

57-63 

T 

ratio  of  specific  heats 

64-70 

X 

blade  length,  see  Figure  38 

Card  3 

1-7 

c 

blade  chord,  see  Figure  38 

8-14 

Awall 

Afclade 

end  wall  and  stiffener  wetted 
area/blade  wetted  area 

Card  4 

1-7 

N 

«  number  of  pressure  ratios  in 

Card  Set  4.  1.  £.  N  ^  20. 
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Figurs  132:  DA  TA  CA  RD  A  R  RANG  EM  ENT  FOR  CASCADE  LOSS  PREDICTIONS 
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Engine  Stability  Margin  Module  (ENGI )  Input 


Figure  133  displays  the  data  card  arrangement  for  the  Engine 
Stability  Margin  Module.  A  description  of  the  card  input 
follows: 


Card 

Column 

Code 

Explanation 

Card  1 

1-4 

ENGI 

control  card  -  contains  the  word 

ENGI 

Card  2 

1-7 

BPR 

bypass  ratio,  2.  £  BPR  ^  12. 

8-14 

TRMIX 

=  1.  indicates  mixed  flow  engine 

=  2.  indicates  non-mixed  flow 
engine 

15-21 

TRAN 

=  -1.  indicates  engine  decel 

=  indicates  steady  state 

operation 

=  +1.  indicates  engine  accel 

22-28 

TTRAN 

time  at  which  engine  accel  or 
decel  ia  initiated,  seconds 

29-35 

POWER 

if  TRAN  =  -1.,  percent  engine 
power  level  et  end  of  engine 
decel  from  100  percent  power 

if  TRAN  =  0 . ,  percent  steady 
state  power  level 

if  TRAN  =  +1.,  percent  engine  power 
level  at  start  of  engine  accel 
to  100  percent  power 

Card  3 

1-7 

N 

number  of  times  and  area  matches 
to  follow  in  Card  Set  4. 

1.  £  N  ^  20. 

Figure  133:  DA  TA  CARD  ARRANGEMENT  FOR  ENGINE  STABILITY  MARGIN  MODULE 


APPENDIX  I 


CHANT'S  THEORY  FOR  THE  ROLLUP  OF  A  JET  IN  CROSSFLOW 

The  complex  potential  for  flow  around  a  circular  cylinder 
is  given  by 

(1) 

where  w  ■  is  the  complex  potential  for  uniform  flow, 

and  the  second  term  is  a  doublet  of  strength  yi/=tTrueBa 
located  at  the  origin  with  its  axis  in  the  -x  direction,  and 
"a"  is  the  radius  of  the  cylinder.  The  complex  variables 
w  and  z  are  given  by 


w  -  ^  il P 


(2) 


?sx+Aj  =  re*e 

(3) 

where  p  and  V  are  the  potential  and  stream  functions.  The 
complex  velocity  field  is  given  by 


*£  =  u-av  (4) 

di 

From  Equations  1  and  2  the  complex  potential  w  is  separated 
into  its  real  and  imaginary  parts  p  and  y>  . 

(5) 


(6) 


The  velocity  components  are  given  by 

U«(l -£)»,«  17) 

+  (8) 
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The  circulation  strength  for  a  segment  of  the  circle  can  be 
replaced  by  a  discrete  vortex  filament,  as  shown  on  the 
following  sketch: 


The  circulation  is  found  by  integrating  the  velocity  vector 
along  a  path  from  point  1  to  point  2. 


Integration  of  Equation  10  gives  the  vortex  strength: 
Pw=  4  ITU*,  asi« 

Referring  to  the  following  sketch: 
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the  velocity  components  induced  at  a  point  (x^,  Yp)  by  a 
vortex  of  strength  r*n  at  the  point  (x^,  y  )  are  given  by 


U 


n 


'£Tr\r 


sine 


rv  yP-y« 

Z7r  (*P-xM)z  +  (y P-y»f 


(12) 


-Ik 

Itrr 


cos© 


pjj.  x?-** 

2tr 


(13) 


The  deformation  of  the  jet  cross  section,  comprised  of  the 
vortex  filaments,  is  found  by  assuming  that  the  filaments 
are  free  to  move  in  space  for  a  finite  time  period.  Their 
direction  of  travel  is  in  the  direction  of  the  velocity  vector 
induced  by  all  the  other  elements.  Thus,  the  lateral  dis¬ 
placement  of  each  vortex  filament  is  found  by: 


Ax  - 

AtZ 

(14) 

n 

Ay  • 

At£v* 

n 

(15) 

where  At  is  the  time  interval  that  is  fixed  by  the  jet 
velocity  and  the  distance  chosen  for  steps  between  cross 
sections,  As. 

The  procedure  is  similar  for  finding  the  deformation  of  a 
cross  section  other  than  circular,  the  difference  being  the 
strength  of  the  vortex  filaments  must  be  determined  from 
some  other  source.  The  three-dimensional  potential  flow 
program  TEA-230  (Ref.  7)  was  used  to  find  the  vortex  strengths 
for  TYPE  =  2.  rectangular  jets.  The  effect  of  rectangle 
aspect  ratio  on  the  vortex  filament  strengths  is  shown  in 
Figure  Al. 

Linear  interpolation  for  vortex  strengths  at  aspect  ratios 
other  than  a/b  -  0.25,  1.0,  and  4.0  was  found  inadequate. 

This  was  not  unexpected  because  linear  superposition  of 
solutions  is  valid  only  if  the  geometry  is  fixed  and  flow 
conditions  vary.  However,  an  asymptotic  equation  for  vortex 
strength  as  a  function  of  a/b  provided  an  excellent  fit  to 
the  data  at  all  intermediate  values. 


r 


C, 
a/ \o  ~ 


+  C3 


(16) 
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VORTEX  STRENGTH. 


0.6 


t 


Figure  A 1:  EFFECT  OF  RECTANGLE  ASPECT  RA  TIO  ON  VORTEX  STRENGTH 
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where  the  constants  C^,  C 2,  and  are  solved  for  by  satisfying 
equation  (16)  at  a/b  -  0.25/  1.0  and  4.0. 


[o.is-c*  ‘  ,-cJ 


(17) 


_  0.2S-4K 
l-K 


(18) 


<=3 


C, 

o.ts-ct 


(19) 


where 


K  ^  0.15 


r,.o-  r 4,,o 
^oi  s  “  r,a 


(20) 


The  flow  direction  relative  to  the  cross  section  has  an  important 
effect  for  rectangular-shaped  sections.  The  effect  of  wind 
direction  is  shown  in  Figure  A2.  Solutions  for  any  wind 
direction  can  be  obtained  by  linearly  combining  solutions  for 
0  and  90  degrees  because  the  geometry  is  fixed. 


r  *  ro„  co zp  -t  P?0.  s\*p 


(21) 


The  vortex  strengths  for  a/b  «  0.25,  1.0,  ar.d  4.0  at  relative 
wind  directions  of  0  and  90  degrees  are  tabulated  in  subroutine 
RCROSS.  Vortex  strengths  for  a  particular  case  are  found  by 
asymptotic  interpolation  (equations  16-20)  to  account  for 
rectangle  aspect  ratio  and  by  Equation  21  to  account  for 
relative  wind  direction. 
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eo 
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Figure  A2:  EFFECT  OF  CROSSWiND  ANGLE  ON  VORTEX  STRENGTH 


APPENDIX  II 


Program  Sample  Cases 


This  appendix  contains  card  input  formats  and  computer  print¬ 
outs  for  sample  cases  computed  by  programs  TEM-356A,  TEM-356B, 
and  TEM-357. 

Jet  Trajectory  and  Spreading  Program  TEM-356A 

A  typical  sample  case  computed  by  TEM-356A  is  described  below. 

TYPE  =1.  Input — This  sample  case  is  for  the  circular  jet 
shown  in  Figure  63.  Program  inputs  are  shown  in  Figure  Bl. 

The  jet's  initial  position  (x,  y,  z)  is  located  at  the  reference 
coordinate  system  origin,  vector  angle  is  135  degrees,  and 
dynamic  pressure  ratio  cj, j/^eo  =  2.78. 

TYPE  =  1.  Printout — Computer  printout  is  given  in  Figure  B2. 
Printout  includes  program  inputs,  {x,  y,  z)  coordinate  out¬ 
put  of  points  on  the  cross  section,  and  summary  data  for  the 
jet  trajectory,  thickness  (DEL),  and  width  (H) .  The  (x,  y,  z) 
coordinates  were  punched  on  cards  for  computer  plotting. 

Reingestion  Prediction  Program  TEM-356B 

Two  sample  cases  were  computed  by  TEM-356B.  Input 
formats  and  computer  printout  for  each  case  are  described 
in  the  following  paragraphs. 

Sample  Case  Number  One  Input — This  sample  case  is  for  an 
early  design  version  o£  the  747  airplane  with  long  duct 
nacelles  and  annular  target  thrust  reversers  as  sketched 
in  Figure  74.  Program  inputs  are  shown  in  Figure  B3. 

Nacelles  are  located  at  30  and  50  percent  span  locations. 

The  inboard  thrust  reverser  plume  is  represented  by  a 

TYPE  «  4.  annular  jet  defined  by  MR  *  10.  arms  vectored 

120  degrees  from  the  3  axis.  Flow  coditions  are  70 

knots  freestream  velocity  (  118.16  ft/sec)  and  40.96. 

Sample  Case  Number  One  Printout — Computer  printout  is  given  in 
Figure  B4.  Printout  includes  inlet  and  jet  input  data  followed 
by  reingestion  diagnostic  predictions.  On  page  five  of  the 
printout,  the  program  predicts  far-field  fountain  reingesticn 
of  the  fourth  jet  arm  into  the  outboard  inlet.  NIP  is  the 
number  of  points  intersecting  an  inlet  streamtube  at  each 
aximuth  angle,  On  page  ten,  crossflow  reingestion  is 

predicted  into  the  outboard  inlet.  Data  summarizing  the 
number  of  points  intersecting  the  inlet  streamtube  are  printed 
on  page  11.  Eight  intersection  points  occur  for  the  third  arm. 
The  (x,  y,  z)  coordinates  of  the  inlet  streamtube  are  printed 
out  for  clock  angles  from  0  degrees  to  360  degrees. 
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Figure  B1:  SAMPLE  CASE  FOR  TYPE  =  1.  CIRCULAR  JET 


JET  I KA  j kC TOR Y  AND  SPREADING  PROGRAM  TEM356A 


DATE  OE  RUN  APR  22,  1972 


fYBt  *  i.  INDICATES  A  CIRCULAR  JfcT  PlUMF 

Ur'c  =  2.  indicates  a  suuake  or  rectangular  jet  plume 

TfPc  =  J •  INDICATES  A  T wO-Dl MENSIUNAL  JET  PLUME 
TYPfc  =  4.  INDICATES  AN  ANNULAR  OR  UMBRELLA  SHAPED  JET 


InPuI  CONDITIONS  FUR  THE  FIRST  JET  FOLLOW  ON  THE  NEXT 


flfrur*  B2:  PROGRAM  PRINTOUT  FOR  TYPE  -  1  SAMPLE  CASE 


PLUME 

PAGf 
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*oun 

TV  PE 
XY20 

TXYZ  * 
QKAT1Q  = 

i 

UINF 

AL  PH  A  = 

PS  I 

> 

D  * 

delsd 

SEND 

I 

I 


» 


0.1E+01 . 

0.0.  0 .0 »  o.c» 

-0.707107E*00,  0.0.  -0.  7071 0 7fc*00, 

0.2  7?778E+0 1 . 

0.5E+02. 

0.0. 

0.0. 

0.2E+01. 

0.2E+00, 


» 


Fi9un  tt:  PROGRAM  PRINTOUT  FOR  TYRE  -  1  SAMPLE  CASE  { Coni ) 


225 


COORDINATE  OUIPUT  Of-  JET 


CROSS  SECTION  NUMBER  1 


POINT 

X 

l 

•  7071 

2  - 

.7011 

3 

.6830 

4 

.6533 

5 

.6124 

6 

.5610 

7 

.5000 

ft 

.4305 

9 

•  3536 

10 

.2706 

11 

.1830 

12 

.0923 

13 

-.0000 

14 

-.0923 

15 

-.1830 

16 

-.2706 

17 

-.3536 

18 

-.4305 

19 

-.5000 

20 

-.5610 

21 

-.6124 

22 

-.6533 

23 

-.6830 

24 

-.7011 

25 

-.7071 

26 

-.7011 

27 

-.6630 

28 

-.6533 

29 

-.6124 

30 

-.5610 

31 

-.5000 

32 

-.4305 

33 

-.3536 

34 

-.2706 

35 

-.1830 

36 

-.0923 

37 

.0000 

38 

.0923 

39 

•  1830 

40 

.2706 

41 

•  3536 

42 

•  4305 

43 

•  5000 

44 

•  5610 

45 

.6124 

46 

.6533 

47 

.6830 

.  68 

.7011 

49 

.7071 

IS  LISTEO  BEL  CM 


V 

l 

0.0000 

-.7071 

-.1305 

-.7011 

-.2588 

-.6830 

-.3827 

-.6533 

-.5000 

-.6124 

— . 6 C88 

-.5610 

-•7C71 

-.5000 

-.7934 

-.4305 

-.8660 

-,,3536 

-.9239 

-r.  2706 

-.9659 

-,1830 

-.9914 

-.0923 

•l.OCOO 

.0000 

-.9914 

.0423 

-.9659 

.1830 

- .9239 

.2706 

-.8660 

.3536 

-.7934 

.4305 

-• 7C7 1 

.5000 

-.6C88 

.5610 

-  •  5COO 

.6124 

-.3827 

.65 33 

-.2588 

.6830 

-.1305 

.7011 

•  0C00 

.7071 

.1305 

.7011 

.2588 

.6830 

.3827 

.6533 

.5000 

.6124 

.6088 

.5610 

.7071 

.5000 

.7934 

.4305 

.8660 

.3536 

.9239 

.2706 

.9659 

.1830 

.9914 

.0923 

1.0000 

-.0000 

.9914 

-.0923 

.9659 

- •  1830 

.9239 

-.2706 

.  8  660 

-.3536 

.7934 

-.4303 

.7071 

-.5000 

.6068 

-.5610 

.5000 

-.6124 

.3827 

-.6533 

.2568 

-.6830 

.1305 

-.7011 

-.0000 

-.7071 

XC 

YC 

1.0000 

0.0000 

.9914 

.1305 

.9659 

.2566 

.92  39 

.3827 

.  8660 

.5000 

.  7934 

.6088 

.7071 

.7071 

•  6  088 

.7934 

.500  0 

.8660 

.1627 

.92  39 

.2  5  66 

.96  59 

.1305 

.9914 

-.0000 

1.0000 

-.1305 

.9914 

-.2588 

.9659 

-.3827 

.92  39 

-• 5000 

.6660 

-.6068 

.7934 

-.7071 

.7071 

-.7934 

.6088 

-.8660 

.5000 

-.9239 

.3827 

-.9659 

.2588 

-.9914 

.1305 

-1.0000 

-.0000 

-.9914 

-.1305 

-.9659 

-.2588 

-.9239 

-.3827 

-.6660 

-.5000 

-.  7934 

-.6088 

-.7071 

-.7071 

-.6088 

-.7934 

-.5000 

- • 8660 

-.3627 

-.9239 

-.2588 

-.9659 

-.1305 

-.9914 

•  0000 

-1.0000 

.1305 

-.9914 

.2588 

-.9659 

.3827 

-.9239 

.5  000 

-.8660 

.6  088 

-.7934 

.7071 

-.7071 

.7934 

-.6088 

•  8660 

-.5000 

.9239 

-.3827 

.9659 

-.2588 

.9914 

-.1305 

I.UOCO 

.0000 

'f&MW  y+CQMH  mWrOUTFO*  TYPE*  1  SAMPLE  CASE  (tint) 


»( 


I 


I 


I 


I 


I 


c 


{. 


POINT 

X 

1 

.4307 

2  . 

•4265 

3 

•  4140 

4 

.3930 

5 

.3635 

6 

.3254 

7 

•  2786 

8 

•  2232 

9 

.1593 

10 

.0874 

11 

•  0081 

12 

-.0777 

13 

-.1688 

14 

-.2637 

1  5 

-.3608 

16 

-.4581 

17 

-.5535 

18 

-.6447 

19 

-.7294 

20 

-.8056 

21 

-.8710 

22 

-.9240 

23 

-.9630 

24 

-.9868 

25 

-.9949 

26 

-.9868 

27 

-.9630 

28 

-.9240 

29 

-.8/10 

30 

-.8056 

31 

-.7294 

32 

-.6447 

33 

-.5535 

34 

-.4581 

35 

36 

-.3608 

-.2637 

37 

-.1688 

38 

-.0777 

39 

•  0081 

40 

.0874 

41 

*1593 

42 

•  2232 

43 

.2786 

44 

•  3254 

45 

.3635 

46 

.3930 

47 

.4  MO 

48 

•  4265 

49 

.4307 

IS  LISTED  BELQn 


V  l 

.0000  -.9860 

-.1122  -.9609 

-•22 36  -.9666 

-.3332  -.9479 

-.4400  -.9189 

-.5429  -.8814 

-.6406  -.8353 

-.7313  -.7808 

-.8136  -.7179 

-.8856  -.6472 

-.9453  -.5692 

-.9909  -.4848 

-1.0206  -.3951 

-I.C329  -.301  7 

- 1 .0264  -.2061 

*1.0003  -.1104 

-.9541  -.0165 

-.8881  .0732 

-.8028  .1366 

-.6997  .2313 

-.5806  .2959 

-.4479  .3461 

-.3047  .3864 

-.1542  .4099 

.0000  .417a 

.1542  .4099 

.3047  .3864 

.4479  .3461 

.5806  .2959 

.6597  .2313 

.8028  .1566 

•8881  .0732 

•9541  -.0165 

1.0003  -.1104 

1.0264  -.2061 

1.0329  -.3017 

1*0206  -.3951 

.9909  -.4848 

.9453  -.5692 

.8856  -.6472 

•8136  -.7179 

.7313  -.7808 

.6406  -.8353 

.5429  -.8814 

.4400  -.9169 

.3332  -.9479 

.2236  -.9686 

*1122  -.9809 

-.0000  -.9850 


XC 

YC 

.63  75 

.0000 

.8316 

.1100 

.8141 

.2191 

.7846 

.3265 

.7433 

.4312 

.6898 

.5320 

.0241 

.62  76 

0463 

.  716b 

.4567 

.797  2 

00  59 

.8677 

•  *446 

.9262 

.1243 

.9709 

-.0035 

1 1.0000 

-. 1 oo6 

1.0120 

-.2730 

1.0057 

-.4095 

.9801 

-.5433 

.9349 

-.6712 

.8701 

-.7901 

.7866 

-• 6*  70 

.6855 

-.9868 

•  5688 

-1.  Ooil 

.4389 

-1.1178 

.2966 

-1.1513 

.1511 

-1.1625 

-.0000 

-1.1513 

-.1511 

-1.1178 

-.2986 

-1.0631 

-.4389 

-.9606 

-.5688 

-.6970 

-.6855 

-.7901 

-.7866 

-•o712 

-.8701 

-.5433 

-.9349 

-.4095 

-.9801 

-.2730 

-1.0057 

-.1368 

-1.0120 

-.0035 

-l .ooo e 

.1243 

-.9709 

.2446 

-.9262 

.3559 

-.8677 

.4567 

-.7972 

•  5463 

-.7166 

.6241 

-.6276 

.6896 

-.5320 

.7433 

-.4312 

•  7646 

-.3265 

.6141 

-.2191 

*8316 

-.1100 

.63  75 

.0000 

*****  mNTOVT FOK7YH.1  SAMPLE  CASE  (Com) 
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COORDINATE  OUTPUT  Of  JET  CROSS  SECTION  NUMBER  SO 

IS  LISTED  BELOW 


POINT 

X 

V 

l 

13.7434 

.0000 

2 

14.0449 

-1.8078 

3 

13.8784 

-3.0268 

4 

14.1111 

-3.6819 

5 

14.137C 

-3.0114 

6 

14.0615 

-2.8328 

7 

14.0125 

-4.0825 

8 

13.9702 

-2.7933 

9 

13.9550 

-1.5421 

10 

14.1246 

-1.8788 

11 

14.2373 

-2.2503 

12 

14.2823 

-3.3561 

13 

14.1982 

-4.2322 

14 

14.0804 

-4.7358 

15 

13.8974 

-3.6122 

16 

13.8763 

-1.7716 

17 

14.0886 

-.7856 

18 

14,3406 

-1.2633 

19 

14.4721 

-1.7672 

20 

14.5562 

-2.4072 

21 

14.575S 

-3.5216 

22 

14.5225 

-4.3263 

23 

14.4387 

-4.3679 

24 

14.3282 

-4.8566 

25 

13.7391 

.0000 

26 

14.3282 

4.8566 

27 

14.4387 

4.3679 

26 

14.5225 

4.3263 

29 

14.5758 

3.5216 

30 

14.5562 

2.4C72 

31 

14.4721 

1.7672 

32 

14.3406 

1.2633 

33 

14.0886 

.7856 

34 

13.8763 

1.7716 

35 

13.8974 

3.6122 

36 

14.0804 

4.7358 

37 

14.1982 

4.2322 

34 

14.2823 

3.3581 

39 

14.2373 

2.2503 

40 

14.1246 

1.8786 

41 

13.9550 

1.5421 

42 

13.9702 

2.7933 

43 

14.0125 

4.0825 

44 

14.0615 
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SUMMARY  Of  JET  CENTERLINE  AND  SPREADING  COEFFICIENT  DATA 


» 


I 


> 


I 


1 


CROSS 


SECT ICN 

S/C 

X/D 

2/0 

XT 

9T 

27 

OEL 

M 

1 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2.0000 

2.0000 

2 

.2000 

-.1410 

.1418 

-.2621 

0.0000 

.2021 

2.0000 

2.0658 

3 

.4000 

-.2797 

.2859 

-.5554 

0.0000 

.5594 

2.0000 

2.1315 

4 

.6000 

-.4131 

.4349 

-.8262 

0.0000 

.8262 

2.0000 

2.1973 

5 

.8000 

-.5369 

.5919 

-1.0739 

0.0000 

1.0739 

2.0000 

2.2631 

6 

l.OOOC 

-.8440 

.  7606 

-1.2801 

0.0000 

1.2001 

2.COOO 

2.3208 

7 

1.2000 

-.7213 

.944  7 

-1.4426 

0.0000 

1.4426 

2.0000 

2.3946 

e 

1.4000 

-.7474 

1. 1421 

-1.4947 

0.0000 

1.4947 

2.0000 

2.4604 

9 

1.6000 

-.7040 

1.3363 

-1.4081 

0.0000 

1.4081 

2.0000 

2.5261 

to 

1.8000 

-.6005 

1 «  5066 

-1.2010 

0.0000 

1.2010 

'2.0000 

2.5919 

it 

2.0000 

-.4607 

1.6492 

-.9213 

0.0000 

.9213 

2.0000 

2.6577 

12 

2.2000 

-.3012 

1.7697 

-.6023 

0.0000 

.6023 

2.0000 

2.7234 

13 

2.4000 

-.130! 

1.873d 

-.2610 

0.0000 

.2610 

2.0000 

2.789:' 

14 

2.6C0C 

•  C471 

1.9680 

•  U941 

0.0000 

-.0941 

2.1073 

3.0622 

IS 

2.800C 

•  2*92 

2 . 04  b  3 

.  4  5  e3 

O.OOOC 

-.4563 

2.4174 

3. 3b44 

16 

3.0000 

.4144 

2.1238 

.62  89 

0.0000 

-.8289 

2.6534 

3.7147 

17 

3.200C 

.6020 

2.1932 

1.2u40 

0.0000 

-1.2040 

2.6949 

4.0529 

18 

3.4000 

.7513 

2.257o 

1.5027 

O.OOOC 

-1.5027 

>. 1418 

4.3986 

19 

3.600C 

.5820 

2. Jldu 

l.So4C 

0.0000 

-1.9640 

3.3939 

4.7514 

20 

J.600C 

1.1738 

2.374r 

2. >4  75 

0.0000 

-2.3475 

3.6508 

5.1112 

21 

4.0000 

1.3664 

2.4285 

2.7i«0 

O.OOOC 

-2.7328 

3.9126 

5.4777 

22 

4.2000 

1.5598 

2.4795 

3.1196 

O.OOOC 

-3.1196 

4.1790 

5.8506 

23 

4.4000 

1.7538 

2.5262 

3. 5075 

O.OOOC 

-3.5075 

4.4499 

6.?29b 

24 

4.6000 

1.5403 

2.5740 

3.6965 

0.0000 

-3.0965 

4.7251 

6.6151 

2S 

4.8000 

2.1432 

2.6195 

4.2664 

O.OOOC 

-4.2064 

5.0045 

7.0063 

26 

5.00C0 

2.3305 

2.662  5 

4.67  71 

0.0000 

-4.6771 

5.2000 

7.4033 

27 

5.20CC 

2.5342 

2.7039 

5.06  84 

0.0000 

-5.0604 

5.5756 

7 .6058 

2b 

5.40CC 

2.7302 

2.7439 

5.4603 

0.0000 

-5.4603 

5.6992 

7. ‘789 

24 

5.60C0 

2.9264 

2.7826 

5.6520 

0.0000 

-5.8528 

5.7680 

8.0753 

30 

5.8000 

3.1220 

2.8201 

6.2457 

0.0000 

-6.2457 

5.0352 

8.1693 

31 

6.0000 

3.3195 

2.8565 

6.6350 

o.onoo 

-6.6390 

5.9009 

8.2612 

32 

6.2000 

3.5163 

2.0919 

7.  0327 

C.0000 

-7.0327 

5.9651 

0.3511 

33 

6.4000 

3.7134 

2.9263 

7.4267 

0.0000 

-7.4267 

6.0279 

8.4391 

34 

6.6000 

3.9105 

2.9550 

7.6211 

0.0000 

-7.0211 

6.0094 

8.5252 

35 

6.8000 

4. 1 C  78 

2.9925 

8.2157 

0.0000 

-0.2157 

6.1497 

6.6096 

36 

7.0000 

4.3053 

3.0243 

0.6106 

0.0000 

-0.6106 

6.2088 

8.6923 

37 

7.2000 

4.5029 

3.0555 

9.0057 

0.0000 

-9.0057 

6.2668 

0.7735 

3b 

7.4000 

4.7005 

3.0059 

9.4010 

0.0000 

-9.4010 

6.3237 

8.0532 

39 

7.600C 

4.0903 

3.1157 

9.7966 

0.0000 

-9.7966 

6.3796 

8.9315 

40 

7.8000 

5.0962 

3.1440 

10.1923 

0.0000 

-10.1923 

6.4345 

9.0004 

41 

8.0000 

5.2541 

3.1734 

1 0.5*02 

0.0000 

-10.5082 

6.4085 

9.0839 

42 

8.2000 

5.4921 

3.2014 

10.9043 

0.0000 

-10.9043 

6.5416 

9.1583 

43 

0.4000 

5.6902 

3.2200 

11.3605 

0.0000 

-11.3005 

6.5930 

9.2314 

44 

0.6000 

5.0804 

3.2557 

11.7769 

0.0000 

-11.77*9 

6.6452 

9.3033 

45 

8.0000 

6.0067 

3.2*22 

12.1733 

0.0000 

-12.1733 

6.6959 

9.3742 

46 

9.0000 

4.2050 

3.3002 

12.5*99 

0.0000 

-12.5699 

6.7457 

9.4440 

47 

9.2000 

6.4033 

3.3337 

12.96*7 

0.0000 

-12.9667 

6.7940 

9.5127 

40 

9.4000 

6.6810 

3.3500 

13.5*36 

0.0000 

•13.3635 

6.0432 

9.5005 

49 

9.6000 

4.0002 

3.3035 

13.76C4 

0.0000 

-13.7604 

6.8909 

9.6473 

50 

9.0000 

7.0707 

3.4079 

14.1575 

0.0000 

-14.1575 

6.9380 

9.7131 

Figun  B2:  PROGRAM  PRINTOUT  FOR  TYPE  -  1  SAMPLE  CASE  (Conclud'd) 
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Figun  B4:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  ONE  (Cont) 
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Figure  B4:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  ONE  (Cont) 
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Figure  B4 :  PROGRAM  PRINTOUT  FOR  TEM456B  SAMPLE  CASE  NUMBER  ONE  { Conti 
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Figure  B4:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  ONE  (Cont) 
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Figure  B4:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  ONE  fCont ) 
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«5jur»fl4:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  TWO  (Concluded) 


This  case  demonstrates  the  need  for  good  arm  definition  for 
TYPE  =  4.  jets.  Because  the  annular  surface  is  represented 
by  discrete  arms,  they  should  be  placed  strategically  pointing 
t  down  and  between  inlets  to  ensure  that  potential  intersection 

points  are  determined  for  far-field  fountain  and  crossflow 
reingestion. 

Sample  Case  Number  Two  Input — This  case  was  designed  to  test 
some  of  the  near-field  fountain  reingesticn  criteria.  Program 
l  inputs  are  given  in  Figure  B5.  Two  inlets  and  three  jets  were 

input. 

Sample  Case  Number  Two  Printout — Computer  printout  is  given  in 
Figure  B6.  Near-field  fountain  reingestion  is  predicted  on 
pages  10  and  11  of  the  printout. 

TR  and  TV  System  Performance  Program 

Thirteen  sample  cases  wore  computed  by  TEM-357.  Input  formats 
and  computer  printouts  for  each  case  are  described  in  the 
following  paragraphs. 

Sample  Case  Number  One  for  ENGI — Program  inputs  and  printout 
are  shewn  in  Figures  B7  and  B8.  Printout  parameters  are 
identified  in  the  following  tabic. 


Parameter 

Explanation 

T 

time 

AMF 

fan  nozzle  area  match 

AMP 

primary  nozzle  area  match 

AMC 

common  nozzle  area  match 

PCF 

thrust  level,  percent 

PCD 

duct  thrust/total  thrust,  percent 

SMF 

fan  stability  margin,  percent 

SMC 

compressor  stability  margin,  percent 

DSMC 

compressor  stability  margin  change  due  to 
engine  acceleration  or  deceleration,  percent 

SMCNET 

net  compressor  stability  margin,  percent 

DPF 

allowable  fan  total  pressure  distortion,  percent 

DPC 

allowable  compressor  total  temperature  distor¬ 
tion,  percent 

DPMAX 

allowable  engine  total  pressure  distortion, 
percent 

DT 

allowable  engine  total  temperature  distortion, 
percent 
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Figun  B6:  PROGRAM  PRINTOUT  FOR  TEM-3S6B  SAMPLE  CASE  NUMBER  TWO 
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Figure  B6:  PROGRAM  PRINTOUT  FOR  TEM-3S6B  SAMPLE  CASE  NUMBER 
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Figure  B6:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  TWO  (Cont) 
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Ffcune  B6:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  TWO  (Cont) 
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Figun  B6:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  TWO  fCont ) 
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Figure  B6:  PROGRAM  PRINTOUT  FOR 
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fflRu/e  B6:  PROGRAM  PRINTOUT  FOR  TEM-356B  SAMPLE  CASE  NUMBER  TWO  (Concluded) 


Figure  B7:  PROGRAM  INPUTS  FOR  ENG /  SAMPLE  CASE  NUMBER  ONE 


sampli  cast  srtditur 
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Figure  B8:  PROGRAM  PR/NOUT  FOR  ENG /  SAMPLE  CASE  FOR  ENGINE  STABILITY  MARGIN  MODULE 


Sample  Case  Number  Two  for  CONI — Program  inputs  and  printout 
are  shown  in  Figures  B9  and  BIO.  Skin  friction,  wetted  area 
(applies  only  if  D^-cl.Q),  underexpansion,  and  bearing  offset 

losses  arc  listed  separately.  The  Reynolds  number  scale  correc¬ 
tion  option  was  used  for  this  case.  Nozzle  Cv  and  CD  were 

corrected  to  a  full  scale  Reynolds  number  Re^  =  10.4x10^. 

Sample  Case  Number  Three  for  ANNU — Program  inputs  and  printout  for 
an  annular  nozzle  are  given  in  Figures  B9  and  Bll.  Skin 
friction,  underexpansion  losses,  and  annular  nozzle  ACV  incre¬ 
ments  relative  to  a  standard  convergent  nozzle  are  listed 
separately. 

Sample  Case  Number  Four  for  EQUI — Program  inputs  and  printout 
appearing  in  Figures  B9  and  B12  are  similar  to  those  for  CONI. 

The  exception  is  a  Cv  penalty  for  increased  wetted  area  due  to 
the  D^e  term. 

Sample  Case  Number  Five  for  TARG — Program  inputs  and  printout 
for  a  clamshell  target  TR  are  given  in  Figures  B13  and  B14. 
Printout  parameters  are  identified  in  the  following  table. 


Parameter 

Explanation 

EBASE 

baseline  corrected  reverser  efficiency 

DELBYD 

reverser  efficiency  increment  due  to  door 

length,  L/D 

DEXBYD 

reverser  efficiency  increment  due  to  setback. 

X/D 

DEHBYD 

reverser  efficiency  increment  due  to  average 

lip  height,  Eh/D 

DESWEP 

reverser  efficiency  increment  due  to  sweep 

angle, ^ 

DEARC 

reverser  efficiency  increment  due  to  arc 

angle,© 

DECONE 

reverser  efficiency  increment  due  to  cone 

angle ,  0( 

DEBEVL 

reverser  efficiency  increment  due  to  bevel 

angle, 0 

ETARC 

corrected  reverser  efficiency,  ^Rc 

ETARG 

static  revei  ency, 

PHI 

airflow  match, 
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Figure  B9:  PROGRAM  INPUTS  FOR  CONI,  ANNU,  &  EQUI 
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fun  810:  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  COMICAL  CRUISE 
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**09912:  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  NONCIRCULAR  NOZZLE 


Fi/un  B13:  PROGRAM  INPUTS  FOR  TARG,  BLOC  SAMPLE  CASE  NUMBER  FIVE.  SIX.  3  SEVEN 


Sample  Case  Number  Six  for  TARG —  Program  inputs  and  printout 
for  an  annular  target  TR  are  given  in  Figures  B13  and  B15. 

Sample  Case  Number  Seven  for  BLOC — Program  inputs  and  printout 
for  a  blocker  external  deflector  TR  are  given  in  Figures  B13 
and  B16.  Printout  parameters  are  identified  in  the  following 
table. 


Parameter 

Explanation 

ERBASE 

baseline  static  reverser  efficiency 

DERBASE 

increment  due  to  blacker  door  angle 

ETARG 

static  reverser  efficiency, 

Kg 

ETARC 

corrected  reverser  efficiency,  77  _ 

1  KC 

PHIBASE 

baseline  airflow  match 

DPHI 

increment  due  to  blocker  door  nngle 

PHI 

airflow  match,  i 

Sample  Case  Number  Eight  for  SING — Program  inputs  and  printout 
for  a  single  bearing  nozzle  are  given  in  Figures  B17  and  B19. 
Theoretical  thrust  components  fcr  this  nozzle  were  given  in 
Figure  29. 

Sample  Case  Number  Nine  ror  THRE — Program  inputs  and  printout 
for  three  bearing  r.ozzle  configuration  number  eight  are  con¬ 
tained  in  Figures  B17  and  B19. 

Sample  Ca.ie  Number  Ten  for  SPHE — Program  inputs  and  printout 
for  a  spherical  eyeball  nozzle  arc  given  in  Figures  B17  and 
B20 . 

Sample  Case  Number  Eleven  for  LOBS--Program  inputs  and  printout  ^ 
are  given  in  Figures  B2]  end  B22  for  a  lobstertail  nozzle. 

Sample  Case  Number  Twelve  for  EXTE — Program  inputs  and  printout 
are  given  in  Figures  B21  and  B23  for  a  hinged  external  deflector. 
Correlations  for  hinged  deflectors  were  developed  during  Part 
IC  and  are  discussed  in  Volume  II  of  this  report.  Special  j 

geometric  inputs  required  are  X/D  ,  D  /D,  and  T  as  defined  in 
the  following  sketch:  r  n 
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Figure  B16:  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  BLOCKER  DEFLECTOR  THRUST  REVERSER 
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Figun  BIB:  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  SINGLE  BEARING  NOZZLE 
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ffour»  572  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  THREE  BEARING  NOZZLE 


SAMPLE  CASE  EO*  SPhEKICAL  EVEttALL  VJZZlE 
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f^purv  020;  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  SPHERICAL  EYEBALL  NOZZLE 


Figure  821:  PROGRAM  INPUTS  FOR  LOBS.  EXTE.  CASC  SAMPLE  CASES  NUMBER  ELEVEN.  TWELVE.  8,  THIRTEEN 
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Figure  B22:  PROGRAM  PRINTOUT  SAMPLE  CASE  FOR  LOBSTERTAIL  NOZZLE 
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Sample  Case  Number  Thirteen  for  CASC— -Program  inputs  and  printout 
are  given  in  Figures  B21  and  B24"  for  a  cascade  lattice  with 
blade  entrance  angle  of  15  degrees  and  solidity  c/s  »  1.5. 

Output  parameters  include  exit  Mach  number  M,,  outlet  flow 
angle  o (2*  total  pressure  loss  Ap^/P^ and  CD. 
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